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Abstract. Here shape space is either the manifold of simple closed smooth 
unparameterized curves in R 2 or is the orbifold of immersions from S 1 to R 2 
modulo the group of diffeomorphisms of S 1 . We investige several Riemannian 
metrics on shape space: L 2 -metrics weighted by expressions in length and 
curvature. These include a scale invariant metric and a Wasserstein type 
metric which is sandwiched between two length-weighted metrics. Sobolev 
metrics of order n on curves are described. Here the horizontal projection of 
a tangent field is given by a pseudo-differential operator. Finally the metric 
induced from the Sobolev metric on the group of diffeomorphisms on M 2 is 
treated. Although the quotient metrics are all given by pseudo-differential 
operators, their inverses are given by convolution with smooth kernels. We 
are able to prove local existence and uniqueness of solution to the geodesic 
equation for both kinds of Sobolev metrics. 

We are interested in all conserved quantities, so the paper starts with the 
Hamiltonian setting and computes conserved momenta and geodesies in gen- 
eral on the space of immersions. For each metric we compute the geodesic 
equation on shape space. In the end we sketch in some examples the differ- 
ences between these metrics. 
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1. Introduction — multiple Riemannian metrics on the space of 

CURVES 

Both from a mathematical and a computer vision point of view, it is of great in- 
terest to understand the space of simple closed curves in the plane. Mathematically, 
this is arguably the simplest infinite-dimensional truly nonlinear space. From a vi- 
sion perspective, one needs to make human judgements like 'such-and-such shapes 
are similar, but such-and-such are not' into precise statements. The common the- 
ory which links these two points of view is the study of the various ways in which 
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the space of simple closed curves can be endowed with a Ricmannian metric. From 
a vision perspective, this converts the idea of similarity of two shapes into a quan- 
titative metric. From a mathematical perspective, a Riemannian metric leads to 
geodesies, curvature and diffusion and, hopefully, to an understanding of the global 
geometry of the space. Much work has been done in this direction recently (see for 
example [10, 12, 14, 15, 23]). The purpose of the present paper is two-fold. On 
the one hand, we want to survey the spectrum of Riemannian metrics which have 
been proposed (omitting, however, the Weil-Peterson metric). On the other hand, 
we want to develop systematically the Hamiltonian approach to analyzing these 
metrics. 

Next, we define the spaces which we will study and introduce the notation we 
will follow throughout this paper. To be precise, by a curve we mean a C°° simple 
closed curve in the plane. The space of these will be denoted B e . We will consider 
two approaches to working with this space. In the first, we use parametrized curves 
and represent B e as the quotient: 

B e S Emb(5 1 ,M 2 )/Diff(S* 1 ) 

of the smooth Frechet manifold of C°° embeddings of S 1 in the plane modulo the 
group of C°° diffeomorphisms of S 1 . In this approach, it is natural to consider all 
possible immersions as well as embeddings, and thus introduce the larger space Bi 
as the quotient of the space of C°° immersions by the group of diffeomorphisms of 
S 1 : 

Emb(5' 1 ,M 2 ) — ► Emb{S\R 2 )/ DifffS 1 ) S B e 

n n n 

Imm(S ,1 ,K 2 ) — ► Imm(S' 1 ,R 2 )/Diff(S' 1 ) ^ B t 

In the second approach, we use the group of diffeomorphisms Diff(R 2 ) of the 
plane, where, more precisely, this is either the group of all diffeomorphisms equal 
to the identity outside a compact set or the group of all diffeomorphisms which 
decrease rapidly to the identity. Let A be the unit circle in the plane. This group 
has two subgroups, the normalizer and the centralizer of A in Diff(R 2 ): 

Diff°(R 2 ,A) c Diff(R 2 ,A) c Diff(R 2 ) 

II II 
{y> | ¥>| A = id A } {<p\<p(A) = £} 

Let i e Emb(5 1 ,R 2 ) be the basepoint i{6) = (sin(0),cos(0)) carrying S 1 to the 
unit circle A. The group Diff (R 2 ) acts on the space Emb(S l , R ) of embeddings by 
composition on the left. The action on the space of embeddings is transitive (e.g., 
choose an isotopy between two embedded circles, transform and extend its velocity 
field into a time-dependent vector field with compact support on R 2 and integrate 
it to a diffeomorphism) . Diff°(R 2 , A) is the subgroup which fixes the base point i. 
Thus we can represent Emb(S' 1 ,R 2 ) as the coset space Diff (R 2 )/ Diff (R 2 , A). 

Furthermore Diff°(R 2 , A) is a normal subgroup of Diff (R 2 , A), and the quotient 
of one by the other is nothing other than Diff (A), the diffeomorphism group of 
the unit circle. So Diff(A) acts on the coset space Diff (R 2 )/ Diff (R 2 , A) with 
quotient the coset space Diff (R 2 )/ Diff (R 2 , A). Finally, under the identification of 
Diff (R 2 )/ Diff °(R 2 , A) with Emb(S' 1 ,R 2 ), this action is the same as the previously 
defined one of Diff on Emb(S' 1 ,R 2 ). This is because if c = ipoi g Emb(S' 1 ,R 2 j, 



RIEMANNIAN METRICS ON SPACES OF CURVES, HAMILTONIAN APPROACH 



3 



and V G Diff (R 2 , A) satisfies ip(i(0)) = i(h(6)), h e Diff^ 1 ), then the action of ip 
carries <p to <p o -0 and hence ctoipoipoi = ipoioh = coh. 

All the spaces and maps we have introduced can be combined in one commuta- 
tive diagram: 

Diff(R 2 ) 

I 

Diff(R 2 )/Diff°(R 2 , A) Emb(S\R 2 ) C Imm(S\R 2 ) 

Diff (R 2 ) / Diff (R 2 , A) B e c B t 

See [12] and [8] for the homotopy type of the spaces Imm(5 1 ,R 2 ) and B t . 

What is the infinitesimal version of this? We will use the notation X(R 2 ) to 
denote the Lie algebra of Diff (R 2 ), i.e., either the space of vector fields on R 2 with 
compact support or the space of rapidly decreasing vector fields. As for any Lie 
group, the tangent bundle TDiff(R 2 ) is the product Diff(R 2 ) x X(R 2 ) by either 
right or left multiplication. We choose right so that a tangent vector to Diff(R 2 ) 
at <p is given by a vector field X representing the infinitesimal curve <p \— > p(x, y) + 
eX(ip(x,y)). 

Fix ip e Diff(R 2 ) and let it map to c = p o i e Emb(5 1 ,R 2 ) and to the curve 
C = Im(c) C R 2 on the three levels of the above diagram. A tangent vector to 
Emb(5 1 , R 2 ) at c is given by a vector field Y to R 2 along the map c, and the vertical 
map of tangent vectors simply takes the vector field X defined on all of R 2 and 
restricts it to the map c, i.e. it takes the values Y(6) = X(c(9)). Note that if c is 
an embedding, a vector field along c is the same as a vector field on its image C. 
A tangent vector to B e at the image curve C is given by a vector field Y along C 
modulo vector fields tangent to C itself. The vertical map on tangent vectors just 
takes the vector field X along c and treats it modulo vector fields tangent to c. 
However, it is convenient to represent a tangent vector to B e or £?, at C not as an 
equivalence class of vector fields along C but by their unique representative which 
is everywhere normal to the curve C. This makes TcBi the space of all normal 
vector fields to C C R 2 . 

In both approaches, we will put a Riemannian metric on the top space, i.e. 
Imm(5 , R ) or Diff(R 2 ), which makes the map to the quotient Bi or to a coset 
space of Diff (R 2 ) into a Riemannian submersion. In general, given a diffcomorphism 
/ : A — > B with a surjective tangent map and a metric G a (h,k) on A, f is a 
submersion if it has the following property: first split the tangent bundle to A 
into the subbundlc TA T tangent to the fibres of / and its perpendicular TA 1 - 
with respect to G (called the horizontal bundle). Then, under the isomorphisms 
df : TAq ^ TBf( a -), the restriction of the A-metric to the horizontal subbundle is 
required to define a metric on TBj, independent of the choice of the point a G f (b) 
in the fiber. In this way we will define Riemannian metrics on all the spaces in our 
diagram above. Submersions have a very nice effect on geodesies: the geodesies on 
the quotient space B are exactly the images of the geodesies on the top space A 
which are perpendicular at one, and hence at all, points to the fibres of the map / 
(or, equivalently, their tangents are in the horizontal subbundle). 
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On Diff(R 2 ), we will consider only right invariant metrics. These are given by 
putting a positive definite inner product G(X, Y) on the vector space of vector fields 
to R 2 , and translating this to the tangent space above each diffeomorphism <p as 
above. That is, the length of the infinitesimal curve ip h- > ip + eX o ip is <J G(X, X). 
Then the map from Diff(M 2 ) to any of its right coset spaces will be a Riemannian 
submersion, hence we get metrics on all these coset spaces. 

A Riemannian metric on Imm(S' 1 ,R 2 ) is just a family of positive definite inner 
products G c (h,k) where c £ Imm(S' 1 ,IR 2 ) and h,k <E C°° (S 1 , K 2 ) represent vector 
fields on R 2 along c. We require that our metrics will be invariant under the action 
of Diff(S' 1 ), hence the map dividing by this action will be a Riemannian submersion. 
Thus we will get Riemannian metrics on Bf. these are given by a family of inner 
products as above such that G c (h, k) = if h is tangent to c, i.e., (h(9),cg(9)) = 
where eg := dgc. 

When dealing with parametrized curves or, more generally, immersions, we will 
use the following terminology. Firstly, the immersion itself is usually denoted by: 

c(9) : S 1 -> R 2 

or, when there is a family of such immersions: 

c(9,t) : S 1 x R 2 . 

The parametrization being usually irrelevant, we work mostly with arclcngth ds, 
arclength derivative D s and the unit tangent vector v to the curve: 

ds = \cg\d9 

D s =d e /\c e \ 

V = Cg/\cg\ 

An important caution is that when you have a family of curves c(9, t), then dg and 
d t commute but D s and d t don't because \cg\ may have a ^-derivative. Rotation 
through 90 degrees will be denoted by: 

The unit normal vector to the image curve is thus 

n = Jv. 

Thus a Riemannian metric on B e or Bi is given by inner products Gc(a, b) where a.n 
and b.n are any two vector fields along C normal to C and a, b <G C°°(C, R). Another 
important piece of notation that we will use concerns directional derivatives of 
functions which depend on several variables. Given a function f(x, y) for instance, 
we will write: 

D( x ,h)f or df(x)(h) as shorthand for dt\of(x + th,y). 

Here the x in the subscript will indicate which variable is changing and the second 
argument h indicates the direction. This applies even if one of the variables is a 
curve C e Bi and h is a normal vector field. 

The simplest inner product on the tangent bundle to Imm(5 1 ,R 2 ) is: 

G° c (h,k) = f (h(9),k(9)) - ds. 
Js 1 
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Since the differential ds is invariant under the action of the group Diff(5' 1 ), the 
map to the quotient Bi is a Riemannian submersion for this metric. A tangent 
vector h to Imm(S' 1 ,R 2 ) is perpendicular to the orbits of Diff(S' 1 ) if and only if 
(h(9),v(9)) = 0, i.e. h is a multiple a.n of the unit normal. This is the same 
subbundle as above, so that, for this metric, the horizontal subspace of the tangent 
space is the natural splitting. Finally, the quotient metric is given by 



where L is a positive definite operator on the vector- valued functions h : S 1 — > R 2 . 
The simplest such L is simply multiplication by some function <f> c (9). However, it 
will turn out that most of the metrics involve L's which are differential or pseudo- 
differential operators. For these, the horizontal subspace is not the natural splitting, 
so the quotient metric on B e and Bi involves restricting G c to different sub-bundles 
and this makes these operators somewhat complicated. In fact, it is not guaranteed 
that the horizontal subspace is spanned by C°° vectors (in the sense that the full 
C°° tangent space is the direct sum of the vertical subspace and the horizontal C°° 
vectors). When dealing with metrics on Diff(R 2 ) and vertical subspaces defined by 
the subgroups above, this does happen. In this case, the horizontal subspace must 
be taken using less smooth vectors. 

In all our cases, L^ 1 will be a simpler operator than L: this is because the 
tangent spaces to B e or Bi are quotients of the tangent spaces to the top spaces 
Diff or Imm where the metrics are most simply defined, whereas the cotangent 
spaces to B e or Bi are subspaces of the cotangent spaces of the space 'above'. The 
dual inner product on the cotangent space is given by the inverse operator and 
in all our cases this will be an integral operator with a simple explicit kernel. A final 
point: we will use a constant A when terms with different physical 'dimensions' are 
being added in the operator L. Then A plays the role of fixing a scale relative to 
which different geometric phenomena can be expected. 

Let us now describe in some detail the contents of this paper and the metrics. 
First, in section 2, we introduce the general Hamiltonian formalism. This is, un- 
fortunately, more technical than the rest of the paper. First we consider general 
Riemannian metrics on the space of immersions which admit Christoffel symbols. 
We express this as the existence of two kinds of gradients. Since the energy function 
is not even defined on the whole cotangent bundle of the tangent bundle we pull 
back to the tangent bundle the canonical symplectic structure on the cotangent 
bundle. Then we determine the Hamiltonian vector field mapping and, as a special 
case, the geodesic equation. We detemine the equivariant moment mapping for sev- 
eral group actions on the space of immersions: the action of the reparametrization 
group Diff (S 1 ), of the motion group of R 2 , and also of the scaling group (if the 
metric is scale invariant). Finally the invariant momentum mapping on the group 
Diff(R 2 ) is described. 




a.b.ds. 



All the metrics we will look at will be of the form: 
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Section 3 is then devoted to applying the Hamiltonian procedure to almost local 
metrics: these are the metrics in which L is multiplication by some function $. Let 
P-c = J s i ds be the length of the curve c and let 

k c {0) = (n(e),D s (v)(0)) = (Jcg,c ee )/\ CB \ 3 

be the curvature of c at c(8). Then for any auxiliary function k), we can define 
a weighted Riemannian metric: 

G*{h,k)= { $(e c ,K c (6))-h{6)k(6)-ds. 
Js 1 

The motivation for introducing weights is simply that, for any 2 curves in Bi, the 
infimum of path lengths in the G° metric for paths joining them is zero, see [12, 13]. 
For all these metrics, the horizontal subspace is again the set of tangent vectors 
a(0)n(0), so the metric on Bi is simply 

Gf(a-n,b-n)= [ $(4, n c {0))a{0)b{0) ■ ds. 
Js 1 

We will determine the geodesic equation, the momenta and the sectional curvature 
for all these metrics. The formula for sectional curvature is rather complicated but 
for special $, it is quite usable. 

We will look at several special cases. The weights 

<&(£,*,) = 1 + Ak 2 

were introduced and studied in [12]. As we shall see, this metric is also closely 
connected to the Wasserstein metric on probability measures (see [1]), if wc assign 
to a curve C the probability measure given by scaled arc length. We show that it is 
sandwiched between the conformal metric G l and G* w where <&w = £~ X + j^^k 2 - 
Weights of the form 

were studied in [10] and independently by [18]. The latter are attractive because 
they give metrics which are conformally equivalent to G°. These metrics are a 
borderline case between really stable metrics on B e and the metric G° for which 
path length shrinks to 0: for them, the infimum of path lengths is positive but at 
least some paths seem to oscillate wildly when their length approaches this infimum. 
Another very interesting case is: 

&(e,K) = e- 3 + A\K\ 2 £- 1 

because this metric is scale- invariant. 

A more standard approach to strengthening G° is to introduce higher derivatives. 
In section 4, we follow the Sobolev approach which puts a metric on Imm(5 1 ,M 2 ) 
by: 

r. n 

G™ m ' n (h,k)= / Y,( D ^ D l k ) ds 

r. n 

= / (Lh,k)ds, where L = ^(-l) 4 L>f 

J S 1 ,;_n 
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However, the formulas we get are substantially simpler and L 1 has an elegant 
expression is we take the equivalent metric: 



We apply the Hamiltonian procedure to this metric. Here the horizontal space 
of all vectors in the tangent space TImm(5' 1 ,K 2 ) which are G lmm ' n -orthogonal 
to the reparametrization orbits, is very different from the natural splitting in §3. 
The decomposition of a vector into horizontal and vertical parts involves pseudo 
differential operators, and thus also the horizontal geodesic equation is an integro- 
differential equation. However, its inverse L _1 is an integral operator whose kernel 
has a simple expression in terms of arc length distance between 2 points on the 
curve and their unit normal vectors. 

For this metric, we work out the geodesic equation and prove that the geodesic 
flow is well posed in the sense that we have local existence and uniqueness of 
solutions in Imm(>S 1 ,R 2 ) and in Bi. Finally we discuss a little bit a scale invariant 
version of the metric G lmm '". For the simplest of these metrics, the scaling invariant 
momentum along a geodesic turns out to be the time derivative of \og(l). At this 
time, we do not know the sectional curvature for this metric. 

In the next section 5, we start with the basic right invariant metrics on Diff(M 2 ) 
which are given by the Sobolev if "-inner product on X(R 2 ). 



These metrics have been extensively studied by Miller, Younes and Trouve and their 
collaborators [5, 14, 15, 20]. Since these metrics are right invariant, all maps to coset 
spaces Diff(K 2 ) — > Diff(R 2 )/ff are submersions. In particular, this metric gives a 
quotient metric on Emb(S' 1 , R 2 ) and B e which we will denote by Gf S ' n (h, k). In 
this case, the inverse L _1 of the operator defining the metric is an integral operator 
with a kernel given by a classical Bessel function applied to the distance in K 2 
between 2 points on the curve. We will derive the geodesic equations: they are all 
in the same family as fluid flow equations. We prove well poscdncss of the geodesic 
equation on Emb(/S' 1 ,R 2 ) and on B e . Although there is a formula of Arnold [4] for 
the sectional curvature of any right-invariant metric on a Lie group, we have not 
computed sectional curvatures for the quotient spaces. 

In the final section 6, we study two examples to make clear the differences be- 
tween the various metrics. The first example is the geodesic formed by the set of 
all concentric circles with fixed center. We will see how this geodesic is complete 
when the metric is reasonably strong, but incomplete in most 'borderline' cases. 
The second example takes a fixed 'cigar-shaped' curve C and compares the unit 
balls in the tangent space TcB e given by the different metrics. 






where L = (1 - AA) n , A = d 2 x + <9 2 . 
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2. The Hamiltonian approach 

In our previous papers, we have derived the geodesic equation in our various 
metrics by setting the first variation of the energy of a path equal to 0. Alternately, 
the geodesic equation is the Hamiltonian flow associated to the first fundamental 
form (i.e. the length-squared function given by the metric on the tangent bundle). 
The Hamiltonian approach also provides a mechanism for converting symmetries of 
the underlying Ricmannian manifold into conserved quantities, the momenta. We 
first need to be quite formal and lay out the basic definitions, esp. distinguishing 
between the tangent and cotangent bundles rather carefully: The former consists of 
smooth vector fields along immersions whereas the latter is comprised of 1-currents 
along immersions. Because of this we work on the tangent bundle and we pull 
back the symplectic form from the contantent bundle to rimm(5' 1 , R 2 ). We use 
the basics of symplectic geometry and momentum mappings on cotangent bundles 
in infinite dimensions, and we explain each step. See [11], section 2, for a detailed 
exposition in similar notation as used here. 

2.1. The setting. Consider as above the smooth Frechet manifold Imm(S' 1 ,R 2 ) 
of all immersions S 1 — > R 2 which is an open subset of C°°(S 1 ,R 2 ). The tangent 
bundle is TImm(S' 1 ,M 2 ) = Inim^ 1 , R 2 ) x C°° (S 1 ,M?) , and the cotangent bundle is 
T* Imm(5' 1 ,K 2 ) = Imm(S ,1 ,R 2 ) x V{S 1 ) 2 where the second factor consists of pairs 
of periodic distributions. 

We consider smooth Ricmannian metrics on Imm(5 1 , R 2 ), i.e., smooth mappings 

G : Imm(S' 1 ,IR 2 ) x C°°(5' 1 ,]R 2 ) x C^S^R 2 ) -» R 

(c, h, k) i — ► G c (h, k), bilinear in h, k 

G c (h,h)>0 for/i^O. 

Each such metric is weak in the sense that G c , viewed as bounded linear mapping 

G c : T c Imm(5 ,1 ,IR 2 ) = C°°(5 1 ,R 2 ) T c * Inim^ 1 , R 2 ) =V(S 1 ) 2 

G : Tlmm(S\M. 2 ) -> T* imm(S' 1 ,R 2 ) 

G(c,h) = (c,G c (h, .)) 

is injective, but can never be surjective. We shall need also its tangent mapping 

TG : T(TImm(S' 1 ,R 2 )) -» T(T* Imm(5' 1 ,M 2 )) 

We write a tangent vector to TImm(5 1 , R 2 ) in the form (c, h;k, £) where (c,h) G 
TImm(S' 1 , R 2 ) is its foot point, k is its vector component in the limits' 1 , Re- 
direction and where i is its component in the C°°(S 1 , Redirection. Then TG is 
given by 

TG(c,h;kJ) = (c,G c (h, .);k,D (c . k) G c (h, .) + G c {£, .)) 

Moreover, if X = (c, h; k, t) then we will write X± = k for its first vector component 
and X 2 = £ for the second vector component. Note that only these smooth functions 
on Imm(S' 1 ,R 2 ) whose derivative lies in the image of G in the cotangent bundle 
have G-gradients. This requirement has only to be satisfied for the first derivative, 
for the higher ones it follows (see [9]). We shall denote by Cg 3 (limits' 1 , R 2 )) the 
space of such smooth functions. 
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We shall always assume that G is invariant under the reparametrization group 
Diff(S' 1 ), hence each such metric induces a Riemann-metric on the quotient space 
^(S 1 , R 2 ) = Imm(S ,1 ,M 2 )/Diff(S* 1 ). 

In the sequel we shall further assume that that the weak Riemannian metric G 
itself admits G-gradients with respect to the variable c in the following sense: 



D c>m G c (h,k) = G c (m,H c (h,k)) = G c (K c (m,h),k) 



where 



H,K : Imm(5 1 ,M 2 ) x G 00 ^ 1 ^ 2 ) x G°°(S\R 2 ) -> G°°(S\R 2 ) 

(c, h, k) 1 ► H c (h, fc), K c (h, k) 

are smooth and bilinear in h, k. 

Note that H and K could be expressed in (abstract) index notation as gij,kg kl and 
9ij,k9 tl - We will check and compute these gradients for several concrete metrics 
below. 



2.2. The fundamental symplectic form on TImm(S' 1 ,R 2 ) induced by a 
weak Riemannian metric. The basis of Hamiltonian theory is the natural 1- 
form on the cotangent bundle T* Imm(5' 1 ,R 2 ) given by: 

9 : T(T* Imm(5 1 ,M 2 )) = imm(S' 1 ,R 2 ) x P(S' 1 ) 2 x G 00 ^ 1 ,^ 2 ) x ^(S 1 ) 2 -» R 

(c, a; h, (3) 1— » (a, h). 

The pullback via the mapping G : TImm(S' 1 , M 2 ) -> T* Imm(5 1 , M 2 ) of the 1-form 
is then: 

(G*e)^ h) (c,h;k,£) = G c (h,k). 

Thus the symplectic form ui = —dG*<d on TImm(5 1 ,R 2 ) can be computed as 
follows, where we use the constant vector fields (c, h) ^ (c, h;k,£): 

W(c,fc)((fel^l),(fe^2)) = -d(G*e)((fei^i),(Ai2,^))|(c,fc) 

= -D {cM) G c (h, k 2 ) - G c (h,k 2 ) + D {cM) G c (h, h) + G c (e 2 , fci) 
(1) = G c (k 2 ,H c (h, h) - K c (k u h)) + G c (£ 2 ,h) - G c (£ u k 2 ) 



2.3. The Hamiltonian vector field mapping. Here we compute the Hamilton- 
ian vcctorficld grad"(/) associated to a smooth function / on the tangent space 
TImm(S\R 2 ), that is / e Cg ) (Imm(S' 1 , M 2 ) x C°°(S' 1 ,M 2 )) assuming that it has 
smooth G-gradicnts in both factors. See [9], section 48. Using the explicit formulas 
in 2.2, we have: 

w (Cih) (grad"(/)( C , h), (k, £)) = uj (c>h) ((grad? (/ )(c, h), grad£(/)(c, h)), (k, I)) = 
= G c (k,H c (h, grad? (/) (c, h))) - G c (K c (grad^ (/) (c,h),h), k) 
+ G c (£, grad^ (/)(c, h)) - G c (grad£(/)(c, h), k) 

On the other hand, by the definition of the cj-gradicnt we have 

w (Cih) (grad"(/)(c, ft), (fc, *)) - df(c, h)(k, i) = D {c . k) f(c, h) + D (M) /(c, h) 
= G c (gradf (/)(c, h), fc) + G c (grad^(/)(c, h), £) 



10 PETER W. MICHOR, DAVID MUMFORD 

and we get the expression of the Hamiltonian vcctorficld: 



gradft/Xc, h) = &ad£(f)(c,h) 

grad£(/)(c, h) =-gradf (/)(c, h)+H c (h, grad^(/)(c, h)) -K c (grad^(/)(c, ft), h) 

Note that for a smooth function / on TImm(S' 1 ,R 2 ) the w-gradient exists if and 
only if both G-gradients exist. 



2.4. The geodesic equation. The geodesic flow is defined by a vector field on 
TImm(S' 1 , R 2 ). One way to define this vector field is as the Hamiltonian vector 
field of the energy function 



E(c,h) = ^G c (h,h), 
The two partial G-gradients are: 

G c {&adg(E)(c,h),e) 
grad$ (E)(c,h) 
G c (gradf(£)(c,/i),fc) 



E : Imm(S' 1 ,M 2 ) x G°°(S\R 2 ) -» 

d 2 E(c,h)(£) = G c (h,£) 
h 

dif?(c,ft)(fc) - \D {Ctk) G c {h,h) 
\G c {k,H c {h,h)) 
\H c {h,h). 



gradf(£)(c,/i) 
Thus the geodesic vector field is 

gr&df (E)(c,h) = h 

grad^Xc, h) = ±H c {h, h) - K c {h, h) 
and the geodesic equation becomes: 
h 



Ct 



ht = \H c {h,h) - K c (h,h) 



or 



ctt = \H c {c t ,c t ) - K c (ct,c t ) 



This is nothing but the usual formula for the geodesic flow using the Christoffel 
symbols expanded out using the first derivatives of the metric tensor. 



2.5. The momentum mapping for a G-isometric group action. We consider 
now a (possibly infinite dimensional regular) Lie group with Lie algebra q with a 
right action g i— > r 9 by isometrics on Imm(S' 1 , M 2 ). If X(Imm(S' 1 , M 2 )) denotes the 
set of vector fields on Imm(S 1 , R 2 ), we can specify this action by the fundamental 
vector field mapping £ : q — > X(Imm(S' 1 , M 2 )), which will be a bounded Lie algebra 
homomorphism. The fundamental vector field (x, X £ $ is the infinitesimal action 
in the sense: 

Cx(c) = 9 t |or cxp(tX) (c). 
We also consider the tangent prolongation of this action on TImm(5 1 , R 2 ) where 
the fundamental vector field is given by 

CI? 1 " 1111 : (c,h) -> (c,h;Cx(c),D (c>h) ({ x )(c) =: C x (c,h)) 

The basic assumption is that the action is by isometries, 

G c (h,k) = ((r 9 )*G) c (h, k)=G rS(c) (T c (r 9 )h,T c {r 9 )k). 
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Differentiating this equation at g = e in the direction X € g we get 

(1) = D {cXx{c)) G c (K k) + G c (C' x (c, ft), k) + G c (h, C' x (c, k)) 

The key to the Hamiltonian approach is to define the group action by Hamiltonian 
flows. To do this, we define the momentum map j : g — > Cg r, (TImm(S' 1 , K 2 ), R) 
by: 

j x (c,h) = G c (Cx(c),h). 

Equivalently, since this map is linear, it is often written as a map 

J : Tlmm(S\K. 2 ) -> fl', (J(c, h),X) = j x (c, ft). 

The main property of the momentum map is that it fits into the following commuta- 
tive diagram and is a homomorphism of Lie algebras: 

H (T Imm) U~ C% (T Imm, R) — ^ X(T 1mm, u) ^ H 1 (T Imm) 






where X(TImm, u>) is the space of vector fields on Tlmm whose flow leaves ui fixed. 
We need to check that: 

Cx(c) = grad^(j x )(c,/i) = &udg{j x ){c,h) 

C x (c, h) = grad^0x)(c, h) = - gradf (j x )(c, h) + H c (h, ( x (c)) - K c (( x (c), h) 

The first equation is obvious. To verify the second equation, we take its inner 
product with some k and use: 

G(fc,gradf (j x )(c, h)) = D (c . k) j x (c, h) = D {c , k) G c (( x (c), h) + G c (C^(c, k), h) 

= G c (k,H c (( x (c),h))+G c (£ x (c,k),h). 

Combining this with (1), the second equation follows. Let us check that it is also a 
homomorphism of Lie algebras using the Poisson bracket: 

{jx ,jv}(c,h) = dj Y (c,h) (grad" (j x ) (c, h) , grad£ (j x ) (c, ft) ) 

= dj Y (c,h)(( x (c),( x (c,h)) 

= D (cXx{c)) G c (CY(c),h) + G c (C Y (cXx(x)),h)+G c (C Y (c),C' x (c,h)) 
= G c (t Y (c,(x(c))-C x (c,(Y(c)),h) by (1) 
- G c ([Cx,Cr](c),ft) = G c (C[x,Y](c),h) = j [x . Y] (c). 

Note also that J is equivariant for the group action, by the following arguments: 
For g in the Lie group let r 9 be the right action on Imm(S' 1 ,M 2 ), then T(r 9 ) 0(^0 
(r 9 ) -1 = CA&ig-^x- Since r 9 is an isometry the mapping T(r 9 ) is a symplectomor- 
phism for ui, thus grad" is equivariant. Thus j x o T(r 9 ) = j\d(g)x phis a possible 
constant which we can rule out since j x (c, h) is linear in ft. 

By Emmy Noether's theorem, along any geodesic t c(i, . ) this momentum 
mapping is constant, thus for any X £ g we have 

(J(c, c t ), X) = j x (c, c t ) = G c (( x (c), c t ) is constant in t. 



We can apply this construction to the following group actions on Imm(S' 1 ,t 2 ). 
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• The smooth right action of the group Diff(S' 1 ) on Inim(S' 1 , R 2 ), given by 
composition from the right: c i— » cop for ip £ Diff(S' 1 ). For X £ X(S 1 ) the 
fundamental vector field is then given by 

CF ff (c) = Cx(c) = 3 t | ( C oFlf) = Ce .X 

where Fl t denotes the flow of X. The reparametrization momentum, for 
any vector field X on S 1 is thus: 

jx(c,h) = G c (c g .X, h). 

Assuming the metric is reparametrization invariant, it follows that on any 
geodesic c{6, t), the expression G c (c$.X,Ct) is constant for all X. 

• The left action of the Euclidean motion group M(2) = M 2 x SO (2) on 
Imm(S' 1 ,M 2 ) given by c ^ e aJ c + B for (B,e aJ ) £ R 2 x 5*0(2). The 
fundamental vector field mapping is 

C(B,o)(c) =aJc + B 

The linear momentum is thus G c {B,h),B £ R 2 and if the metric is trans- 
lation invariant, G c (B,c t ) will be constant along geodesies. The angular 
momentum is similarly G c (Jc,h) and if the metric is rotation invariant, 
then G c (Jc,c t ) will be constant along geodesies. 

• The action of the scaling group of R given by c ^ e r c, with fundamental 
vector field ( a (c) — a.c. If the metric is scale invariant, then the scaling 
momentum G c (c, c t ) will also be invariant along geodesies. 

2.6. Metrics and momenta on the group of diffeomorphisms. Very similar 
things happen when we consider metrics on the group Diff(R 2 ). As above, the 
tangent space to Diff (R 2 ) at the identity is the vector space of vector fields X(R 2 ) 
on R 2 and we can identify TDiff(R 2 ) with the product Diff(R 2 ) x X(R 2 ) using 
right multiplication in the group to identify the tangent at a point p with that 
at the identity. The definition of this product decomposition means that right 
multiplication by ip carries (ip, X) to (<poip, X). As usual, suppose that conjugation 
p tpopoip^ 1 has the derivative at the identity given by the linear operator Ad^, 
on the Lie algebra X(R 2 ). It is easy to calculate the explicit formula for Ad: 

Ad^A) = (Dip ■ X) o ip -1 . 

Then left multiplication by ip on Diff(M 2 ) x X(R 2 ) is given by (<p,X) i— > (ip o 
(p,Ad^(X)). We now want to carry over the ideas of 2.5 replacing the space 
Imm(5 1 ,R 2 ) by Diff(R 2 ) and the group action there by the right action of Diff (R 2 ) 
on itself. The Lie algebra g is therefore X(M 2 ) and the fundamental vector field 
Cx(c) is now the vector field with value 

(x(<f) = dt\o(<P H P cxp(tA) o p- 1 ) = Ad v (X) 

at the point p. We now assume we have a positive definite inner product G(X, Y) on 
the Lie algebra X(M 2 ) and that we use right translation to extend it to a Riemannian 
metric on the full group Diff (M 2 ). This metric being, by definition, invariant under 
the right group action, we have the setting for momentum. The theory of the last 
section tells us to define the momentum mapping by: 

jx(<P,Y) = G(Oc(<p),Y). 
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Noether's theorem tells us that if (p(t) is a geodesic in Diff (ffi 2 ) for this metric, then 
this momentum will be constant along the lift of this geodesic to the tangent space. 
The lift of <p(t), in the product decomposition of the tangent space is the curve: 

hence the theorem tells us that: 

G(Ad v ( t )(X),d t ((fi) o <£ _1 (t)) = constant 
for all X. If we further assume that Ad has an adjoint with respect to G: 

G(Ad v (X),Y) = G(X,KA* V (Y)) 
then this invariance of momentum simplifies to: 



Ad*( t ) (dt(ip) ° <p 1 (t)) = constant 



This is a very strong invariance and it encodes an integrated form of the geodesic 
equations for the group. 

3. Geodesic equations and conserved momenta for almost local 

riemannian metrics 

3.1. The general almost local metric G*. We have introduced above the $- 
metrics: 

Gf(h,k):= f <t>(£ ClKc (d))(h(0),k(0))ds. 
Js 1 

Since £(c) is an integral operator the integrand is not a local operator, but the 
nonlocality is very mild. We call it almost local. The metric G* is invariant under 
the reparametization group Diff (S 1 ) and under the Euclidean motion group. Note 
(see [12], 2.2) that 

D {c , h) e c = [ { -^ld9= [ (D s (h),v)ds 

= — I (h, D s (v))ds = — / n(c)(h,n}ds 
Js 1 Js 1 

(Jh 9 ,c e e) (Jc g ,h gg ) (h ,c e } 

D( c ,h)Kc = — j— To — + — j— 3k(c 

Icel 3 \c e \ 2 

= (D 2 s (h),n)-2 K (D s (h),v). 
We compute the G*-gradients of c i— > G*(/i, k): 

D {c>m) G*(h, k) = J K).D (c>m) e c .{h, k) + d 2 <S>(£, n).D {c . m) K c .(h, k) 

+ $(e,K).{h,k){.D s (m),v)) ds 

= — I K c (m,n)ds- / di$(£, n)(h,k)ds 
Js 1 Js 1 

+ J (d 2 $(l,K)({D 2 s (m),n) -2 K {D s (m), «)) + *(4/s)<D a (m),v))<M}d8 
= J $(e,K)(m, K \j d 1 $(l,K)(h,k}ds)n + D 2 s (d 2 <S>{£, K )(h,k)n) 
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+ 2D s (d2$(£,K)n{h,k)v) - D s ($(e, n)(h, k)v H \ds 



According to 2.1 we should rewrite this as 

D (c . m) Gt(h, k) = G*(K?(m, h), k) = G*(m, H?(h, k)), 
where the two G*-gradients if* and if* of c i— » Gf(h, k) are given by: 



Kf(m,h) = —(^J K c (m,n)d, 



+ (<^(m), n) - 2«<I> 8 (m) ) «)) ft + (D s (m), v)h 

H*(h, k) = ^J-^ {-{ Kc J fc > ds ) n + D * ( 92$ ^' fc H + 

+ 2i? s k)k(H, £>) - D s (&(t, n)(h, fc)u) 

By substitution into the general formula of 2.4, this gives the geodesic equation 
for G*, but in a form which doesn't seem very revealing, hence we omit it. Below 
we shall give the equation for the special case of horizontal geodesies, i.e. geodesies 
in Bi. 

3.2. Conserved momenta for G*. According to 2.5 the momentum mappings 
for the reparametrization, translation and rotation group actions are conserved 
along any geodesic 1 1— > c(t, . ): 



$>(l c , k c ){v, c t )\ce\ 2 G X(5 1 ) reparametrization momentum 

/ $(^ c , K c )c t ds eR 2 linear momentum 

Js 1 

/ §(£ c , k c )(Jc, c t )ds e M angular momentum 
Js 1 



Note that setting the reparametrization momentum to and doing symplectic 
reduction there amounts exactly to investigating the quotient space Bi(S 1 ,M 2 ) = 
Imm(5 1 , R 2 )/ Diff (S 1 ) and using horizontal geodesies for doing so; a horizontal 
geodesic is one for which (v, c t ) = 0; or equivalently it is G*-normal to the Diff (5* 1 )- 
orbits. If it is normal at one time it is normal forever (since the reparametrization 
momentum is conserved). This was the approach taken in [12]. 

3.3. Horizontality for G*. The tangent vectors to the Diff(S' 1 ) orbit through c 
are T c (c o Diff {S 1 )) = {X.c : X e G°°(S' 1 ,M)}. Thus the bundle of horizontal 
vectors is 

N c = {h e C^iS 1 ^ 2 ) : (h,v) = 0} 

= {a.n e G°°(5' 1 ,IR 2 ) : a E G°°(5' 1 ,]R)} 

A tangent vector h e T c Imm(5 1 , R 2 ) = G°° (S 1 , R 2 ) has an orthonormal decompo- 
sition 

h = h T + h ± G T c (c oDiff + (5 1 ))eAA c where 
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h T = (h, v)v G T c (c O Diff+(5 1 )), 
h 1 - = (h, n)n e W c , 

into smooth tangential and normal components, independent of the choice of <&(t, n) . 
For the following result the proof given in [12], 2.5 works without any change: 

Lemma. For any smooth path c in Imm(S , M. ) there exists a smooth path tp in 
Diff(S' 1 ) with <p(0, . ) = Idgi depending smoothly on c such that the path e given 
by e(t,6) = c(t,tp(t,6)) is horizontal: e t -Leg. □ 

Consider a path t i— > c(-,t) in the manifold Imm(S' 1 ,R 2 ). It projects to a path 
7r o c in Bi(S 1 , R 2 ) whose energy is called the horizontal energy of c: 



For a horizontal path this is just the usual energy. As in [12], 3.12 we can express 
E hol (c) as an integral over the graph S of c, the immersed surface Scl 3 param- 
eterized by (t,6) i > (t,c(t,0)), in terms of the surface area dfis = \®t x &e\d6dt 
and the unit normal ns = (rig, ng, n|) of S: 



Here the final expression is only in terms of the surface S and its fibration over the 
time axis, and is valid for any path c. This anisotropic area functional has to be 
minimized in order to prove that geodesies exists between arbitrary curves (of the 
same degree) in B l (S 1 ,M 2 ). 

3.4. The horizontal geodesic equation. Let c(8, t) be a horizontal geodesic for 
the metric G*. Then c t (9,t) = a(6,t).n(6,t). Denote the integral of a function 
over the curve with respect to arclength by a bar. Then the geodesic equation for 
horizontal geodesies is: 



This comes immediately from the formulas for H and K in the metric G when 
you substitute m = h = k = a.n and consider only the n-part. We obtain in this 
case: 






$ • (K, n) = -(«a)A$.a + d 2 <S>.D 2 s (a).a + <9 2 $.kV - $«a 2 



$ • (H, n) = -{di<$>a 2 ).n + £> 2 (<9 2 $.a 2 ) + <9 2 $.K 2 a 2 - $«a 2 . 
and the geodesic formula follows by substitution. 
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3.5. Curvature on B i j(S 1 ,« 2 ) for G*. We compute the curvature of B i (S 1 ,M 2 ) 
in the general almost local metric G*. We proceed as in [12], 2.4.3. We use the 
following chart near G e B i (6' 1 ,R 2 ). Let c e Imm / (5 1 ,R 2 ) be parametrized by 
arclength with 7r(c) — C of length L, with unit normal n c . We assume that the 
parameter 9 runs in the scaled circle S£ below. 

i, : C°°(Sl (-e, e)) - Im m/ (5i, M 2 ), Q(c) := ^(G°°(Si, (-e, s))) 

^(/)(0) = c(0) + f(6)n c (6) = c(0) + f(6)ic'(6), 

it o ^ : G°°(Si, (-e.e)) - B^^.R 2 ), 

where £ is so small that ^>(/) is an embedding for each /. We have (see [12], 2.4.3) 

v>(/)' = c' + fid + fid' = (i - fn c y + fid 

W)" = d' + fid + 2 fid' + fid" = -(2/'« c + K)c' + (« c + /" - fn 2 c )id 



1 U\-fK c )id-fd\ 



^(l-/«c) 2 +r 

T/V-Zi = h.id G G°°(S' 1 ,M 2 ) = T W) Imm/^.R 2 ) 

Mi-/«c) ft/ 7 



(1 - M) 2 + /' 2 (l-/«c) 2 + / 

(Tfip.h) 1 - = Ml~/«c) g 
v /(l-/« c ) 2 + .f 2 

««/) = ((1 _ fe) l + f2) 3/ 2 ^(/)'W)"> 

= « c + (/" + /«?) + (/ 2 «? + §/' 2 « c + //'4 + 2//"« c ) + o(/ 3 ) 

W/)) = / I W)l d» = / (1 - 2/« c + J 2 * 2 + /' 2 ) 1/2 d9 

= (l - /«c + ^- + 0(/ 3 )) d6 = L-JT c + \f + 0(f) 

where we use the shorthand 5 = J gl d# = J g i ds. Let G* denote also 

the induced metric on Bij(Sj j , M 2 ). Since 7r is a Riemannian submersion, for / G 
C°°(Sl, (-e,e)) and /i,fce G°°(S'i,]R) we have 



((tt o ip)*G*)f(h, k) = Gt Wf)) (T f (n o tfh, T f (n o 
= G* (/) ((T / V./i) ± ,(T / V.fc) ± ) = / *W(/)W(/)) 



hk(l- fn c 



v /(l-/« c ) 2 + /' 2 

We have to compute second derivatives in / of this. For that we expand the main 
contributing expressions in / to order 2: 



(1 - /«) 2 (1 - 2/k + /V + .f' 2 )- 1 / 2 = 1 - /« - if + 0(f) 
k) = k c ) + K C )(* - L) + 9 2 $(L, k c )(k - K C ) 
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+ d 1 d 2 $(L, Kc )(l- L)(K- K C ) 

+ d_ML^ {l _ Lf + dJHL^ {K Kc)2 + 0(3) 

We simplify notation as k — k c , $ = k c ), ((iroip)*G®) / = G* etc. and expand 
the metric: 



Gj(h, k) = ^ hk - di*.fn + &$.(/" + /k 2 ) - 



+ + 9 2 <i>.(/ 2 « 3 + y 2 K + ff K > + 2ffK) 

+ d^.fnj^ - &$./«.(/" + fn 2 ) - $.±/' 2 ) d9 + 0(f) 

Note that Gq(/i, fc) = J gl hk$d9. We differentiate the metric and compute the 
Christoffel symbol at the center / = 

-2G A (T {h, k), I) = -dG A (0)(l)(h, k) + dG A (0)(h)(k, I) + dG A (0)(k)(l, h) 
= J (-d^M.M-diQ.hM.l + diQ.hk J lKd8 1 - d 2 <S>" .hkl 

- 2d 2 &.h'kl - 2d 2 &.hk'l - 2d 2 <5>.h'k'l + d 2 <$>.hkln 2 - ^.hklnj d8 



Thus 



1 



T (h, k) = — [d-L$.(hK.k + h.kn) - ndi^.hk 

+ d 2 $".hk + 2d 2 &.tik + 2d 2 &hk' + 2d 2 $.h'k' 
- d 2 ®.hkn 2 + $.Wck] 



Letting h — k — f t = a, this leads to the geodesic equation from 3.4. For the 
sectional curvature we use the following formula which is valid in a chart: 

2R f (m,h,m,h) = 2Gf(R f (m,h)m,h) = 

= -2d 2 G A {f)(m,h)(h,m) + d 2 G A (.f)(m,m)(h,h) + d 2 G A (,f)(h,h)(m,m) 
- 2G A (T(h, m),T(m, h)) + 2G A {T(m, m),T(h, h)) 

The sectional curvature at the two-dimensional subspace Pf(m,h) of the tangent 
space which is spanned by m and h is then given by: 

, % % G*(R(m,h)m,h) 
kf(P f (m,h)) = — ' 



|m|| 2 ||/i|| 2 -G*(m,h) 2 ' 

We compute this directly for / = 0, using the expansion up to order 2 of G A (h, k) 
and the Christoffels. We let W(6 1 , 6 2 ) = h{e 1 )m(e 2 )-h(e 2 )m(e- L ) so that its second 
derivative d 2 W{6 u 6»i) = W 2 (0i,0i) = h{e{)m'{9{) - /i/(0i)ra(0i) is the Wronskian 
of h and m. Then we have our final result for the main expression in the horizontal 
sectional curvature, where we use J — J sl , g — J sl gds, and $i = <9i$ etc. Also 
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recall that the base curve is parametrized by arc-length. 



i?o (m, h, m, h) = G* (Ro(m, h)m, h) = 

$ $, M - 2f<K) 2 - (§ 2 n) 2 \ 

«.*2 - | + 2 2 ^ (0l)W 2 (01, ^l) 2 d*i 

$22(^1) , 



+ 



+ 



/ (^^-^§§-^) (^1)^2(^1,^1) / w(0 1 ,9 2 ) K (e 2 )de 2 de 1 

- $12) (0l)W 22 (01, ft ) / 2 )k(0 2 ) <#2 <#1 



k(0 2 )k(0 3 Wi, 2 )W(0i, 3 ) d0 2 Mi d9 3 



3.6. Special case: the metric G A . If we choose $(£ c ,k c ) = 1 + Ak 2 then we 
obtain the metric used in [12], given by 

G*(h,k)= [ (l + AK c (6) 2 )(h(9),k(9))ds. 
Js 1 

As shown in our earlier paper, \/I is Lipschitz in this metric and the metric domi- 
nates the Frechet metric. 

The horizontal geodesic equation for the G^-metric reduces to 

-\n c a 2 + A (a 2 (-D 2 ( Kc ) + |/t|) - 4D 8 (K c )aD s (a) - 2n c D s {a) 2 ) 



a t 



1 + 



as found in [12] ,4.2. Along a geodesic t c(t, . ) we have the following conserved 
quantities: 

(1 + AkI){v, c t )\cg\ 2 G X(5' 1 ) reparametrization momentum 

/ (1 + An 2 )c t ds € R 2 linear momentum 

Js 1 

/ (1 + Ak 2 )( Jc, c t )ds € R angular momentum 

Js 1 

For k) = 1 + An 2 we have <9i$ = 0, d 2 $ = 2Ak, df $ = 2A, and the general 
curvature formula in 3.5 for the horizontal curvature specializes to the formula in 
[12], 4.6.4: 

n *, , ^ f/ C 1 - Ak 2 ) 2 -4A 2 kk" + 8A 2 k ,2 j „ 2 . 2 \ ja 
i?*(m, h,m,h)= J ( - i 2 (1 + ^ 2 ) W 2 2 + AW 2 2 2 J 
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3.7. Special case: the conformal metrics. We put <&(f (c), k(c)) — &(£(c)) and 
obtain the metric proposed by Menucci and Yezzi and, for <f> linear, independently 
by Shah [18]: 

G?(M) = $(4) / (h,k)ds = $(l c )G° c (h,k). 
Js 1 

All these metrics are conformally equivalent to the basic L 2 -metric G°. As they 
show, the infimum of path lengths in this metric is positive so long as $ satifies an 
inequality <&(t) > C.£ for some C > 0. This follows, as in [12], 3.4, by the inequality 
on area swept out by the curves in a horizontal path c t = a.n: 

J \a\M<(J a 2 .d s y%/ 2 < (_y 1/2 . (G * (a7a)) i/ 2 

a + + ^ ( ^(t,-) \ 1/2 , x . G*-path length 
Area swept out < max ——7 — — ■ (G -path length < — . 

- * \*(tc(t,)J Vc 

The horizontal geodesic equation reduces to: 

« 2 W /i/f 2j \ f /* 
a t = —a — 1 t=t 1 I a .as ] k — \ / n.a.ds 



2 $ 



If we change variables and write b(s,t) = $(l(t)).a(s,t), then this equation simpli- 
fies to: 



2$ I $ 



Along a geodesic 1 1— » c(t, . ) we have the following conserved quantities: 
$(^ c )(u, c t )\c'{6)\ 2 £ X(S 1 ) reparametrization momentum 

$>(£ c ) / Ctds e K 2 linear momentum 

^(^c) / ( Jc, c t )rfs G R angular momentum 

For the conformal metrics, sectional curvature has been computed by Shah [18] 
using the method of local charts from [12]. We specialize formula 3.5 to the case 
that &(1,k) = is independent of k. Then 82$ = 0. We also assume that h,m 
are orthonormal so that $/i 2 = $m 2 = 1 and $/im = 0. Then the the sectional 
curvature at the two-dimensional subspace Po{m, h) of the tangent space which is 
spanned by m and h is then given by: 

G* (i?o(m, h)m, h) 



k (P {m, h)) = 



||m||W-Gj(m,fc)2 



2$ ( + ] 4$3 
which is the same as the equation (11) in [18]. Note that the first line is positive 
while the last line is negative. The first term is the curvature term for the H - 
metric. The key point about this formula is how many positive terms it has. This 
makes it very hard to get smooth geodesies in this metric. For example, in the case 
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where = c.£, the analysis of Shah [18] proves that the infimum of G* path 
length between two embedded curves C and D is exactly the area of the symmetric 
difference of their interiors: Area(Int(C)A Int(D)), but that this length is realized 
by a smooth path if and only if C and D can be connected by 'grassfire', i.e. a 
family in which the length \c t {0,t) \ = 1. 



3.8. Special case: the smooth scale invariant metric G SI . Choosing the 
function k) = fr z + A^j- we obtain the metric: 

-2 , 



Gf(M) = (±+A^)(h,k)ds. 



The beauty of this metric is that (a) it is scale invariant and (b) log(^) is Lipschitz, 
hence the infimum of path lengths is always positive. Scale invariance is clear: 
changing c, h, k to A • c, A • h, A • k changes £ to A • £ and k to k/A so the A's in G SI 
cancel out. To sec the second fact, take a horizontal path ct = a- n, < t < 1, and 
abbreviate the lengths of the curves in this path, £ c (t,-)i to £(t). Then we have: 

d\og£(t) _ 1 



/ K c{t,-)(9) ' a(9,t)ds, hence 
Js 1 



dt £{t) 
8\og£{t) ffK 2 a 2 ds\ 1/2 fjl-ds^ 1/2 



£{t) ) v m 

< -j= (G SI (a, a)) 1/2 , hence 



dt 



| log(^(l)) - log(f(0))| < Si-path length/ VA 

Thus in a path whose length in this metric is K, the lengths of the individual curves 
can increase or decrease at most by a factor e K l^. Now use the same argument 
as above to control the area swept out by such a path: 

J \a\ds < (^J a 2 d^j ■ ^ 1 ■ ds^j 

< (£ 3 G SI {a,a)) 1/2 ■ l 1/2 = £ 2 ■ G SI {a,a) 1/2 , hence 

Area-swept-out < e K/VJ £(0) 2 ■ K 

which verifies the second fact. We can readily calculate the geodesic equation for 
horizontal geodesies in this metric as another special case of the equation for G*: 

at = l + AWf ( ^ + A{iK) ^ - M " D ^y ^£ 2 nD s (a) 2 

q 4/2 \ 

- 4:A£ 2 D s (n)aD s (a) + (3 + A(£k) 2 ) 

where the "overline" stands now for the average of a function over the curve, i.e. 
/ ■ ■ - ds/£. Since this metric is scale invariant, there are now four conserved quan- 
tities, instead of three: 

<&(f, k)(v, ct)\c'(9)\ 2 e X(5 1 ) reparametrization momentum 

/ <S>(£,K,)c t ds e M 2 linear momentum 

Js 1 
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I 

L 



k)(Jc, ct)ds £ K angular momentum 

s 1 



<&(£, k)(c, c t )ds € R scaling momentum 

s 1 

It would be very interesting to compute and compare geodesies in these special 
metrics. 

3.9. The Wasserstein metric and a related G*-metric. The Wasserstein met- 
ric (also known as the Mongc-Kantorovich metric) is a metric between probability 
measures on a common metric space, see [1], and [2] for more details. It has been 
studied for many years globally and is defined as follows: let p and v be 2 proba- 
bility measures on a metric space (X, d). Consider all measures p on X x X whose 
marginals under the 2 projections are p and v. Then: 



d w&ss (p,v) = inf // d(x,y)dp(x,y). 

p--Pi,*(P)=H,P2,*(P)=v J J XxX 

It was discovered only recently by Benamou and Brcnicr [6] that, if X = R™, this is, 
in fact, path length for a Riemannian metric on the space of probability measures 
V . In their theory, the tangent space at p to the space of probability measures and 
the infinitesimal metric are defined by: 

T^.-p = | vector fields h = V/ completed in the norm J \h\ 2 dp}^ 

where the tangent h to a family p(t) is defined by the identity: 

^ + div(M = o. 

In our case, we want to assign to an immersion c the scaled arc length measure 
p c = ds/i. This maps Bi to V. The claim is that the pull-back of the Wasserstein 
metric by this map is intermediate between G l and G* w , where 

$w(£,k) =£~ 1 + T5^ 2 - 

This is not hard to work out. 

(1) Because we are mod-ing out by vector fields of norm 0, the vector field h 
is defined only along the curve c and its norm is i~ x . j \\h\\ 2 ds. 

(2) If we split h = av + bn, then the condition that h = V/ means that 
J a.ds = and the norm is l~ x . J (a 2 + b 2 )ds. 

(3) But moving c infinitcsimally by h, scaled arc length parametrization of c 
must still be scaled arc length. Let c(-,t) = c + th. Then this means 
\cg\ t — enst. \c$\ at t = 0. Since \cg\ t = (c t e, c$)/\c$\, this condition is the 
same as (D s (av + bn), v) = enst., or D s a — bn c = enst.. 

(4) Combining the last 2 conditions on b, we get a formula for a in terms of b, 
namely a — K * (bn c ), where we convolve with respect to arc length using 
the kernel K(x) = sign(a;)/2 - x/l, -I < x < i. 

(5) Finally, since \K * /|(o;) < \K\.\f \ = yJTJVl\!\ for all /, it follows that 

r 1 . J b 2 d s < r 1 . J (a 2 + b 2 )d s < r 1 . J (b 2 + ^--b 2 )ds 
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which sandwiches the Wasserstein norm between G e and G*"' for $vk 
I- 1 . {1 + {Ik) 2 / 12). 



4. Immersion-Sobolev metrics on Imm(5 ,R ) and on Bi 

4.1. The G lmm, ™-metric. We note first that the differential operator D s = ^ is 
anti self-adjoint for the metric G°, i.e., for all h, k £ C°° (S 1 M 2 ) we have 

J (D a (h),k)da = J (h,-D a (k))da 

We can define a Sobolev-type weak Riemannian metric 1 on Imm(S' 1 ,R 2 ) which is 
invariant under the action of Diff (S 1 ) by: 

(1) G™(M) = f ((h,k)+A.(Dl l h,D^k)).ds 

Js 1 

= / (L n (h),k)ds where 
Js 1 

(2) L n (h) or L n Jh) = (I + (-l) n A.D 2n )(h) 



The interesting special case n = 1 and A — > oo has been studied by Trouve and 
Younes in [20, 23] and by Mio, Srivastava and Joshi in [16, 17]. In this case, the 
metric reduces to: 

G immT,oo (M) ^ f {Ds{h)jDs ( k]) . ds 

Js 1 

which ignores translations, i.e. it is a metric on Imm(S' 1 ,M 2 ) modulo translations. 
Now identify M 2 with C, so that this space embeds as follows: 

Imm(S' 1 ,M 2 )/transl. ^ C loo (5 1 ,C) 
c i — ► Cg. 

Then Trouve and Younes use the new shape space coordinates Z(6) = ce(6) and 
Mio ct al use the coordinates <&(#) = \og(cg(9)) - with complex square roots and logs. 
Both of these unfortunately require the introduction of a discontinuity, but this will 
drop out when you minimize path length with respect to reparametrizations. The 
wonderful fact about Z{9) is that in a family Z(t, 8), we find: 

Z t = so \ Zt \ 2 d9 = \ J ^r\ce\d6 = \ J \D s {c t )\ 2 ds 

so the metric becomes a constant metric on the vector space of functions Z. With 
$, one has / |$ t | 2 ds = / \D s (c t )\ 2 ds, which is simple but not quite so nice. One 
can expect a very explicit representation of the space of curves in this metric. 

Returning to the general case, for each fixed c of length £, the differential op- 
erator L n c is simply the constant coefficient ordinary differential operator / 



^There are other choices for the higher order terms, e.g. summing all the intermediate deriva- 
tives with or without binomial coefficients. These metrics are all equivalent and the one we use 
leads to the simplest equiations. 
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/+(-l)"A/( 2n ) on the s-line modulo l.Z. Thus its Green's function is a linear com- 
bination of the exponentials exp(A.x), where A are the roots of 1 + (— f )"AA 2 ™ = 0. 
A simple verification gives its Green's function (which we will not use below): 

A 2 " = (-l)«+i/A 

This means that the dual metric G™ m '™ = (G™" 1 '™) -1 on the smooth cotangent 
space C 00 (5 ,1 ,IR 2 ) = G° C (T C Imm(5' 1 , K 2 )) C T* ImmfS 1 ,! 2 ) = V^S 1 ) 2 is given by 
the integral operator L" 1 which is convolution by F n with respect to arc length s: 

G~(h,k) = ff F n ( Sl - s 2 ).(h(si),k(s2)).ds 1 .ds 2 . 

JJs 1 xS 1 



4.2. Geodesies in the G lmm, ™-metric. Differentiating the operator D s = j^dg 

with respect to c in the direction m we get — ^r^^-dg, or —(D s m,v)D s . Thus 
differentiating the big operator L„ jC with respect to c in the direction m, we get: 

2n-l 

(3) D [c . m) L n Ah) = (-1)" +1 A ]T Dl{D s (m),v)D 2 s n -i(h) 

j=o 

Thus we have 



D (Ctm) GT m ' n (h,k) 



2n-l 



= A. f (-1)" +1 ]T (D 3 s (D s m,v)D 2 s n ^ (h),k)ds + f (L n (h),k)(D s m,v)ds 

„ 2n—l „ 

= A.J Y(-l) n+ 3 +1 ((D s m,v}D 2n -3(h),Diky s + J (h,k)(D s m,v}ds 

= j sl ( m ' A ^(-l) n+j D s ((D 2n -ih,Dik)v^ -D s ((h,k)v)y s 

According to 2.1 we should rewrite this as 

D {c . m) G~ l (h,k) = G~(K?(m,h),k) = G™ m ' n {m,H:(h,k)), 

and thus we find the two versions K n and H n of the G n -gradient of c G™ m ' n (/i, k) 
are given by: 

(4) K2(m,h) = L- 1 ((-ir+ 1 A. Di(D s m,v)D 2n -i(h) + (D s m,v)h) 
and by 

(2n-l 
A. J2(-ir+W s ((D 2n -ih,Dik)v) - D s ((h,k)v)) 

3=1 ' 

(2n-l 2i 

3 = 1 3=2 
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2n-l -. 

+ {D 2 s n - ] h,Dik)n c n- (D s h,k)v - (h,D s k)v - (h,k) Kc n) 

=L~ l ^ - (L„(fc), L> s fc)w - (D a /i, L n (k))v - (h, k)n{c)n 

2n-l 

(5) +A.J2(-l) n+i (Ds n - i h,Dik)K(c)n) 
0=1 J 
since D s (y) = n{c)n. By 2.4 the geodesic equation for the metric G n is 



c M = i^ c "( Ct ,c t )-^™(c t ,c t ). 



We expand it to get: 



(6) 



L n {ctt) = -(L n (ct),D a (ct))v- 



\c t \ 2 K{c) 



n - (D a (c t ),v)c t 



A 



2n-l 



+ -.J2(-ir+i{D^c t ,Dic t ) K (c)7 



3 = 1 



+ H) n AE Di((D s (c t ),v)D 2 s n -i(c t )) 

3 = 1 



From (3) we see that 



2n-l 



(L„(ct)) t - L„(c„) = dL n (c)(c t )(c t ) = (-1)" +1 A ^ D^Ct),^ 2 ""^*). 

3=0 

so that a more compact form of the geodesic equation of the metric G n is: 



(7) 



(L n (c t ))t = -(L n (c t ),D s (ct))v - 



\Ct\ 2 K{c) 



n - (D s (c t ),v)L n c t 



A 



2n-l 



+ J- E (-V n+j (D 2n -jc t ,Dic t )K(c)n 



i=i 



For n = this agrees with [12], 4.1.2. 



4.3. Existence of geodesies. 



Theorem. Let n > 1. For eac/i fc > 2n+l i/ie geodesic equation 4-2 (6) /las unique 
local solutions in the Sobolev space of H k -immersions. The solutions depend C°° 
on t and on the initial conditions c(0, . ) and Ct(0, . ). The domain of existence 
(int) is uniform in k and thus this also holds in Imm(5 1 ,K 2 ). 



Proof. We consider the geodesic equation as the flow equation of a smooth (C°°) 
vector field on the # 2 -open set U k x H k (S 1 ,M 2 ) in the Sobolev space H k (S 1 ,M 2 ) x 
H k (S 1 ,R 2 ) where U k = {c e H k : \c e \ > 0} C H k is # 2 -opcn. To see that this 
works we will use the following facts: By the Sobolev inequality we have a bounded 
linear embedding if^S 1 , R 2 ) C C m (S\lR 2 ) if k > m + §. The Sobolev space 
iJ fe (5' 1 ,R) is a Banach algebra under pointwise multiplication if k > |. For any 
fixed smooth mapping / the mapping u i— ► / o u is smooth H k — > if fc > 0. 
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The mapping (c, u) i— > L niC u is smooth £/ x i? fe — » jj k - 2n an d is a biboundcd 
linear isomorphism i? fe — » H k ~ 2n for fixed c. This can be seen as follows (see 
4.5 below): It is true if c is parametrized by arclength (look at it in the space 
of Fourier coefficients). The index is invariant under continuous deformations of 
elliptic operators of fixed degree, so the index of L n c is zero in general. But L n c 
is self-adjoint positive, so it is injective with vanishing index, thus surjective. By 
the open mapping theorem it is then bibounded. Moreover (c,w) <— > L~ c (w) is 
smooth U k x H k ~ 2n — > H k (by the inverse function theorem on Banach spaces). 
The mapping (c, /) ^ DJ = jKdef is smooth H k x H m D U x H m -» H™- 1 

I c | 

for k > m, and is linear in /. Let us write D c f — D s f just for the remainder of 
this proof to stress the dependence on c. We have v = D c c and n = JD c c. The 
mapping c i— ► n{c) is smooth on the i? 2 -open set {c : \cg\ > 0} C H k into H k ~ 2 . 
Keeping all this in mind we now write the geodesic equation as follows: 

c t = u=: Xi(c,u) 

u t = - (L n , c (u), D c (u))D c (c) - |Ctl 2 K(c) JD c (c) - (D c (u),D c c)u 

+ t E (-l) n+3 (^^^^)«(c)^(c) 

2«-l 

+ (-1)"A 2 ^ c ((D c (u),D c (c))Dt n -i(u)) ) 
=: X 2 (c,u) 

Now a term by term investigation of this shows that the expression in the brackets 
is smooth U k x H k — > H k ~ 2n since fc - 2n > I > i. The operator L~* then takes 
it smoothly back to i? fe . So the vector field X = (Xi, X2) is smooth on U k x i7 fe . 
Thus the flow Fl exists on H k and is smooth in t and the initial conditions for 
fixed k. 

Now we consider smooth initial conditions c = c(0, . ) and u = c t (0, . ) = 
u(0, . ) in C°°(S' 1 , M 2 ). Suppose the trajectory Fl t (c , Uo) 01 through these intial 
conditions in H k maximally exists for t e (— a/j, bk), and the trajectory Fl* +1 (c , uo) 
in iJ' c+1 maximally exists for t G (— afc + i,6fc + i) with fe^+i < 6^. By uniqueness 
we have Fl k+1 (c ,u ) = Fl*(co,uo) for t € (— afe+i,6fe+i). We now apply 9e to the 
equation Ut — Xi(c, u) = L~ c ( . . . ), note that the commutator [de, L~ c ] is a pseudo 
differential operator of order — 2n again, and write w = dgu. We obtain wt = dgUt = 
Ln]c®8 (•••) + [®0 1 £«.c] (•••)• I n tnc term d$ ( . . . ) we consider now only the terms 
dg n+1 u and rename them dg n w. Then we get an equation wt(t,8) = X2(t,w(t,6)) 
which is inhomogeneous bounded linear in w € H k with coefficients bounded linear 
operators on H k which are C°° functions of c, u e if fe . These we already know on 
the intervall (— d/s,6fe). This equation therefore has a solution w(t, . ) for all t for 
which the coefficients exists, thus for all t G (ak,bk)- The limit lim t yb k+1 w(t, . ) 
exists in H k and by continuity it equals D c {u) in H k at t = bk+i- Thus the H k+1 - 
flow was not maximal and can be continued. So (— dk+i, bk+i) = (a>k,bk)- We can 
iterate this and conclude that the flow of X exists in (~} m>k H m = C°° . □ 
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4.4. The conserved momenta of G lmm '™. According to 2.5 the following mo- 
menta are preserved along any geodesic t \— > c(t, . ): 



{co,L ntC (ct))\co(6)\ elfS 1 ) 


reparametrization momentum 


/ L n , c (c t )ds = c t ds eR 2 


linear momentum 






/ (Jc,L n . c (c t ))ds e R 

JS 1 


angular momentum 



4.5. Horizontality for G imm ' n . ft e T c Imm(S ,1 ,K 2 ) is G™ m ' n -orthogonal to the 
Diff(S' 1 )-orbit through c if and only if 

O = G j : ~(ft,Cx(c)) = Gj ; ~(ft,c .X)= / X.(L niC (h),ce)da 

Js 1 

for all X e X(5' 1 ). So the G imm <"-normal bundle is given by 

^ = {^r(S,M 2 ) : (i„,c(ft),^) = 0}. 

The G lmm,rl -orthonormal projection T c Imm — » TV™, denoted by ft i— > ft- 1 = ft- 1 ' 6 
and the complementary projection ft ft T e T c (c o Diff(S' 1 )) are determined as 
follows: 

ft T = X(h).v where (L„ ;C (ft), w) = (L ntC (X(h).v),v) 
Thus we are led to consider the linear differential operators associated to L n .c 
Lj,Lj : G°°(S' 1 ) — » G°°(5' 1 ), 

ic(/) = <in,c(/.«),«> = <in,c(/.n),n) ) 

ic (/) - (£n,c(/.«),n) = -<L n , c (/.n),«). 

The operator is of order 2n and also unbounded, self-adjoint and positive on 
L 2 (S 1 ,\c e \d6) since 

/ Lj(f)gds= / (L n . c (fv),v)gds 
Js 1 Js 1 

= / (fv,L n . c (gv))ds = / fLj(g)ds, 
Js 1 Js 1 

f Lj(f)fds= f (fv,L n Jfv))ds > if / ^ 0. 

In particular, Lj is injective. L;f , on the other hand is of order 2n — 1 and a similar 
argument shows it is skew-adjoint. For example, if n = 1, then one finds that: 

Lj = -A.D 2 S + (1 + A.k 2 ).I 

L^r = -2A.k.D s - A.D s (k).I 

Lemma. The operator Lj : G°°(5 1 ) -> G°°(5 1 ) is invertible. 

Proof. This is because its index vanishes, by the following argument: The index 
is invariant under continuous deformations of elliptic operators of degree 2n. The 
operator 

Lj(f) = (-^ n r4^ d e n (f) + lower order tcrms 
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is nomotopic to (1 + (— l) n dg n )(f) and thus has the same index which is zero since 
the operator 1 + (— l) n dg n is invertiblc. This can be seen by expanding in Fourier 
series where the latter operator is given by (/(m)) i— > ((1 + m 2n )f(m)), a linear 
isomorphism of the space of rapidly decreasing sequences. Since Bj is injective, it 
is also surjective. □ 

To go back and forth between the 'natural' horizontal space of vector fields a.n 
and the G lmm, ™-horizontal vector fields {h | (Lh, v) = 0}, we only need to use these 
operators and the inverse of L T . Thus, given a, we want to find b and / such that 
L(an + bv) = fn, so that an + bv is _ff"-horizontal. But this implies that 

Zr^a) = -(L(an),v) = (L(bv),v) = L T (b). 

Thus if we define the operator C c : C 00 ^ 1 ) -» C°°(S' 1 ) by 

C c := (Lj)-'oL^, 

we get a pseudo-differential operator of order -1 (which is an integral operator), so 
that a.n + C(a).v is always i? lmm '™-horizontal. In particular, the restriction of the 
metric G lmm ' n to horizontal vector fields hi — ai.n + bi.v can be computed like this: 

G~(/n,/i 2 )= / (Lh 1 ,h 2 ).ds 
Js 1 

{L(a\.n + b\.v), n).a2-ds 

= f (L T (oi) + L- L (6i)).o 2 .ds 
Js 1 

= / (L T + L^- oC)a 1 .a 2 .ds. 

Thus the metric restricted to horizontal vector fields is given by the pseudo differ- 
ential operator L rod = L T + L 1 - o (L T ) _1 o L^. On the quotient space B i7 if we 
identify its tangent space at C with the space of normal vector fields a.n, then: 



G l c mm '"(ai, 02) = / (L T + L 1 ' o [L 1 )- 1 o L ± )ai ■ a 2 ■ ds 
Jc 



Now, although this operator may be hard to analyze, its inverse, the metric on 
the cotangent space to Bi, is simple. The tangent space to Bi at a curve C is 
canonically the quotient of that of Imm(S' 1 , R 2 ) at a parametrization c of C, modulo 
the subspace of multiples of v. Hence the cotangent space to Bi at C injects into 
that of Imm(S' 1 ,R 2 ) at c with image the linear functionals that vanish on v. In 
terms of the dual basis v, n, these are multiples of n. On the smooth cotangent 
space C°°(5' 1 ,]R 2 ) S G° C (T C Imm(5 1 , R 2 )) C T* Imm(5 1 ,I 2 ) S ^(S* 1 ) 2 the dual 
metric is given by convolution with the elementary kernel K n which is a simple 
sum of exponentials. Thus we need only restrict this kernel to multiples a{s).h c (s) 
to obtain the dual metric on Bi. The result is that: 

G"(ai,a 2 )= // K n (si - S2).(n c (s 1 ),n c (s2)).a 1 (si).a 2 (s2).dsids2. 

JJs 1 xS 1 
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4.6. Horizontal geodesies. The normal bundle Af c mentioned in 4.5 is well de- 
fined and is a smooth vector subbundlc of the tangent bundle. But Imm(S' 1 ,R 2 ) — > 
B^S 1 , M?) = Imm / Diff (S 1 ) is not a principal bundle and thus there are no princi- 
pal connections, but we can prove the main consequence, the existence of horizontal 
paths, directly: 

Proposition. For any smooth path c in Imm(5' 1 , M 2 ) there exists a smooth path tp 
in Diff(S' 1 ) with ip(0, ■ ) = Id s i depending smoothly on c such that the path e given 
by e(t,6) = c(t,tp(t,6)) is horizontal: (L n _ e {et),eg) = 0. 

Proof. Writing D c instead of D s we note that D coip (fotp) = ^Q^^ = (D c (f))op 
for tp g Diff + (5' 1 ). So we have L n , colp (f o tp) = (L n>c f) o tp. 

Let us write e — co <p for e(t, 9) — c(t, tp(t, 9)), etc. We look for ip as the integral 
curve of a time dependent vector field 9) on S 1 , given by ip t — £ o <p. We want 
the following expression to vanish: 

(L n ,co v (d t (c O p)),dg(co tp)) = (L n ^ coip (c t Otp+{cgOtp) tp t ), (eg O tp) tp ) 

= (L n A c t) ol P + L n , c ( c e-0 °<P,ce° <p)<Pe 

= ({{L n , c {Ct), Cg) + (L niC (^.Cg),Cg)) o tp) tpg. 

Using the time dependent vector field £ = — j^(Lj)~ 1 ((L n _ c (c t ), v)) and its flow 
tp achieves this. □ 

If we write 

c t = na + vb = (n, v^j (^j 

then we can expand the condition for horizontality as follows: 
D s (c t ) = (D s a + K(c)b)n + (D s b - k(c)o)v. 



K(c t ) = c t + (-irA(n,v)^ K 



D 2 s -k 2 D s K + K D s \ n fa 
-D s k-kD s D 2 -k 2 J U 



= c t + (-l) n A(n,v) 
so that horizontality becomes 

o = { L nAct ), V ) = M + ( -ir ( o, i) (J;-_ K * Ds % + - K K ? s Y (; 

We may specialize the general geodesic equation to horizontal paths and then 
take the v and n parts of the geodesic equation. For a horizontal path we may 
write L n ^ c (c t ) = an for a(t, 9) = (L niC (c t ),n) . The v part of the equation turns out 
to vanish identically and then n part gives us (because n t is a multiple of v): 



I 1 2 ( \ ( \ 

~a t = _ (D sCt ,v)a+ ^ ]T (-l) n+J (D 2 s n - 3 cuDic t ) 
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Note that applying 4.3 with horizontal initial vectors gives us local existence and 
uniqueness for solutions of this horizontal geodesic equation. 



4.7. A Lipschitz bound for arclength in G lmm ". We apply the inequality of 
Cauchy-Schwarz to the derivative of the length function £{c) — J \cg\d8 along a 
path 1 1 ^ c(t, . ): 

(cte, eg) 

au = 

'Si 



d t £(c) = d£(c)(c t ) 



s 1 \ c e\ 



[ (D s (ct),v) 
Js 1 



ds 



< 



< 



(<4lM| G i, 



< ^/I(cjC(A,n)\\c t \\ G n, 
d t £(c) C(A,n) 

< : I Ct G n - 



2^£{c) ~ 2 
Thus we get 

ww,.))- ^( C (o, . ))i < 1 \a ty /m\ dt < ^4f± jf 

_ C{A,n) 



•Ct || G™ dt 



-L G n( C ). 



Taking the infinimum of this over all paths t i— > c(t, . ) from cq to ci we see that 
for n > 1 we have the Lipschitz estimate: 

Iv^R- V^MI < ^dist^ m ( Cl ,c ) 
Since we have Lq°«( c ) ^ ^G™( c ) with equality for horizontal curves we also have: 



Ifn>l, \V^-VeiCo)\<ldi S t^(C\,C ) 



4.8. Scale invariant immersion Sobolev metrics. Let us mention in passing 
that we may use the length of the curve to modify the immersion Sobolev metric 
so that it becomes scale invariant: 

G™ m ' scal ' n (M) = / U(c)- 3 (h,k) +£(c) 2n - 3 A(D^(h),D^(k))) ds 
Js 1 

= I ((£(c)- 3 + (-l) n £(c) 2n - 3 ADl n )h,k) ds 
Js 1 

This metric can easily be analyzed using the methods described above. In particular 
we note that the geodesic equation on Imm(5' 1 ,]R 2 ) for this metric is built in a 
similar way than that for G lmm '" and that the existence theorem in 4.3 holds for 
it. Note the conserved momenta along a geodesic t c(t, . ) are: 



£(c 



^3/ c t ds + (-l) n £(c) 2n - 3 A [ D 2n (c t )ds 
c ) Js 1 Js 1 



£{cf 



c t ds € R 2 



linear momentum 



£(c 



^ / (Jc,ct)ds + (~l) n £(c) 2n - 3 A [ (Jc,D 2n (c t ))ds angular momentum 
c ) Js 1 Js 1 
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-J-To I (c, c t ) ds + (-l) n £(c) 2n - 3 A I (c,D 2n (c t )) ds scaling momentum 

«( c ) Js 1 Js 1 

As in the work of Trouve and Younes [20, 23], we may consider the following variant. 

G™ m ' scal ' n '°°(/i,fc) = lim 1 / ((£(c)- 3 + (-l) n £(c) 2n - 3 AD 2 s n )h,k)ds 

A^co A J s i 

= (-i) n e(c) 2n - 3 [ (D 2n h, k) ds 

It is degenerate with kernel the constant tangent vectors. The interesting fact is 
that the scaling momentum for (7 lmm > sca i4,°° j s given by 



—±-r f (c,D 2 s (c t ))ds = d t log£(c). 



5. SOBOLEV METRICS ON Diff(M 2 ) AND ON ITS QUOTIENTS 

5.1. The metric on Diff(K 2 ). We consider the regular Lie group Diff(M 2 ) which 
is either the group Diff c (R 2 ) of all diffcomorphisms with compact supports of R 2 or 
the group Diffs(R 2 ) of all diffeomorphisms which decrease rapidly to the identity. 
The Lie algebra is £(R 2 ), by which we denote either the Lie algebra X C (R 2 ) of vector 
fields with compact support or the Lie algebra Xs(R 2 ) of rapidly decreasing vector 
fields, with the negative of the usual Lie bracket. For any n > 0, we equip Diff(R 2 ) 
with the right invariant weak Riemannian metric G Dlff ' n given by the Sobolev H n - 
inncr product on X C (R 2 ). 



i-\-j<n 

= V (-A) i+J —— — — — / (d 2 ld 2 JX,Y)dx 

n^T^„ i\j\(n-i-j)\ J R 2 



0<i,j<n 
i+j<n 

= [ (LX,Y)dx where 
L = L A , n = (1 - AA) n , A = d 2 xl +d 2 2 . 

(We will write out the full subscript of L only where it helps clarify the meaning.) 
The completion of X C (R 2 ) is the Sobolev space H n {R 2 ) 2 . With the usual L 2 - 
inner product we can identify the dual of H n (R 2 ) 2 with _ff~"(R 2 ) 2 (in the space 
of tempered distributions). Note that the operator L : H n (R 2 ) 2 -> H^ n (R 2 ) 2 is a 
bounded linear operator. On L 2 (R 2 ) the operator L is unbounded selfadjoint and 
positive. In terms of Fourier transform we have LA n u(£,) = (I + A\£\ 2 ) n u. Let 
F A . 

n in the space of tempered distributions 5'(R 2 ) be the fundamental solution (or 
Green's function: note that we use the letter 'F' for 'fundamental' because 'G' has 
been used as the metric) of La,u satisfying L,A, n {FA,n) = <5o which is given by 
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The functions Fa,u are given by the classical modified Bessel functions K T (in the 
notation, e.g., of Abramowitz and Stegun [7] or of Matlab) by the formula: 

FAAX) = 2«7r(„ - iy.A • (^) 

and it satisfies (L^ 1 u)(x) = J R2 F(x — y)u(y) dy for each tempered distribution u. 
The function Fa, u is C" _1 except that Fa,i has a log-pole at zero. At infinity, 
Fa,u{x) is asymptotically a constant times x n ~ 3 / 2 e~ x : these facts plus much much 
more can be found in [7]. 

5.2. Strong conservation of momentum and 'EPDiff '. What is the form of 
the conservation of momentum for a geodesic ip(t) in this metric, that is to say, a 
flow x i — ► ip(x, t) on M 2 ? We need to work out Ad^, first. Using the definition, we 
see: 

/ (LX,Ad* v (Y)) := f (LAd v (X),Y)= [ {(dtp.X) o tp~ 1 ,LY) 
Jr 2 ' Jr 2 Jr 2 

= [ det(d<p)(d<p.X,LY o <p) = [ (X, dct(d^).d^ T .(LY o v )) 
Jr 2 Jr 2 

hence: 

Ad* (F) = L- 1 (dct(d(p).dip T .(LYo<p)) . 

Now the conservation of momentum for geodesies tp(t) of right invariant metrics on 
groups says that: 

L- 1 (det(d V )(t).d V (tY. (l^o^o^ 

is independent of t. This can be put in a much more transparent form. First, L 
doesn't depend on t, so we cross out the outer L _1 . Now let v(t) — ^oip^ 1 e X(M 2 ) 
be the tangent vector to the geodesic. Let u(t) = Lv(t), so that: 

det{d(p){t).d(p{tf ,{u{t) o tp(t)) 

is independent of t. We should not think of u{t) as a vector field on M 2 : this is 
because we want (u, v) to make invariant sense in any coordinates whatsoever. This 
means we should think of u as expanding to the differential form: 

Lu(t) = (ui.dx 1 + U2-dx 2 ) ® fl 

where n — dx 1 A dx 2 , the area form. But then: 

ip(t)*(uj(t)) = (dip t .(uoif(t)),dx) ®det(dtp)(t).n 

so conservation of momentum says simply: 



(p(t)*ui(t) is independent of t 



This motivates calling u(t) the momentum of the geodesic flow. As we mentioned 
above, conservation of momentum for a Riemannian metric on a group is very 
strong and is an integrated form of the geodesic equation. To see this, we need only 
take the differential form of this conservation law. v(t) is the infinitesimal flow, so 
the infinitesimal form of the conservation is: 

d 



-w{t) + C v(t) {w{t)) =0 
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where C v u\ is the Lie derivative. We can expand this term by term: 

4,(t)(ui) = £V-|^ 
j 

dv l 



C v{t) {dx l ) = dv l = ^J- dxJ 



C v(t) (uj{t)) ■■ 



^2 ( v j .-^r.dx i +u j .-^-7.dx i ) + divv. ^mdx 1 I <8> A*- 



The resulting differential equation for geodesies has been named EPDiff: 



V = 


dip _j 




u = L(v) 


dui 
dt 




dui 
dxi 


■ dv j \ 

+ v? . — — 7 + divv.Ui = 0. 
dx % ) 



Note that this is a special case of the general equation of Arnold: d t u — — ad(u)*u 
for geodesies on any Lie group in any right (or left) invariant metric. The name 
'EPDiff' was coined by Holm and Marsden and stands for 'Euler-Poincare', although 
it takes a leap of faith to see it in the reference they give to Poincare. 



5.3. The quotient metric on Emb(S' 1 ,R 2 ). We now consider the quotient map- 
ping Diff(IR 2 ) — ► Emb(S' 1 ,]R 2 ) given by ip ip o i as in the section 1. Since this 
identifies Emb(5 1 ,R 2 ) with a right coset space of Diff(K 2 ), and since the metric 
G2 iS is right invariant, we can put a quotient metric on Emb(5 1 ,M 2 ) for which this 
map is a Riemannian submersion. Our next step is to identify this metric. Let 
ip e Diff(M 2 ) and let c = tp o i g Emb(5 1 , K 2 ). The fibre of this map through tp is 
the coset 

tp. Diff°(S' 1 ,R 2 ) = {V | V ° c = c}.ip. 
whose tangent space is (the right translate by tp of) the vector space of vector fields 
X g X(M 2 ) with IoceO. This is the vertical subspace. Thus the horizontal 
subspacc is 



jy J (LY,X)dx = 7 ifXoc=oj 



If we want Y g X(M 2 ) then the horizontal subspace is 0. But we can also search for 
Y in a bigger space of vector fields on R 2 . What we need is that LY = c*(p(6).ds), 
where p is a function from S 1 to M 2 and ds is arc-length measure supported on C. 
To make c*(p(6).ds) pair with smooth vector fields X(M 2 ) in a coordinate invariant 
way, we should interpret the values of p as 1-forms. Solving for Y, we have: 



Y(x) 



Js 1 



c{9)).p{9)ds 



(where, to make Y a vector field, the values of p are now interpreted as vectors, using 
the standard metric on W 2 to convert 1-forms to vectors). Because F is not C°°, 
we have a case here where the horizontal subspace is not given by C°° vector fields. 
However, we can still identify the set of vector fields in this horizontal subspace 
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which map bijectively to the C°° tangent space to Emb(5' 1 ,]R 2 ) at c. Mapped to 
T c Emb(5' 1 ,M 2 ), the above Y goes to: 

(Y o c)(6) = f F(c(9) - c(6» 1 )).p(0 1 ).|c'(^ 1 )|d6» 1 
Js 1 

(1) =:(F c *p)(6) where 



F c (6 1 ,6 2 ) = F(c(6 1 )-c(6 2 )) = ±- f 

^ Jr 



_ / r i(c(0i)-c(e 2 ),e 1 dr 



Note that here, convolution on S 1 uses the metric L 2 (S' 1 , \c'(9\d9) and it defines a 
self-adjoint operator for this Hilbcrt space. Moreover, it is covariant with respect 
to change in parametrization: 

F co <p * (/ oip) = (F c * /) o tp. 

What are the properties of the kernel F c l From the properties of F, we see that 
F c is C™ -1 kernel (except for log poles at the diagonal when n = 1). It is also a 
pseudo-differential operator of order — 2n + 1 on 5 1 . To see that let us assume for 
the moment that each function of 9 is a periodic function on R. Then 

c(0i) - c(9 2 ) = f c e (9 2 + - 9 2 ))dt.(9 x - 9 2 ) =: c(9 1} 9 2 ){9 X - 9 2 ) 
Jo 



F (9, 9o) = — f e l{ - 01 - 02){ ~ c{01 > e2) ^ - df 

1 f M,-e^( f \5(0i,0 2 )\ 



-2 



di]2^J drji 



=: J_ f e ^- e ^F c (9 u 92,m)dvi 

where we changed variables as 771 = (c(6\, 6*2), £) and 772 = ( Jc(6\, 62), 0- So we see 
that F c {9\,92) is an elliptic pseudo differential operator kernel of degree — 2n + 1 
(the loss comes from integrating with respect to 772). The symbol F c is real and 
positive, so the operator p 1— > F c *p is self-adjoint and positive. Thus it is injective, 
and by an index argument similar to the one in 4.5 it is invcrtible. The inverse 
operator to the integral operator F c is a pseudo-differential operator L c of order 
2n — 1 given by the distribution kernel L c (9,9i) which satisfies 

L c *F c *f = F c *L c *f= f 

(2) L cov *(hoip) = ((L c *h)oip) for all tp e Diff+(S' 1 ) 

If we write h — Y o c, then we want to express the horizontal lift Y in terms of 
h and write Yh for it. The set of all these Yh spans the horizontal subspace which 
maps isomorphically to T c Emb(S' 1 , M 2 ). Now: 

h = Y o c = (F * (c* (p.ds))) o c = F c * p. 

Therefore, using the inverse operator, we get p = L c * h and: 

Y h = F* (c*(p.ds)) = F*(c*((L c *h).ds)) or 

Y h (x)= [ F{x-c{9))[ L c {9,9 1 )h{9 1 )\c'{9 1 )\d9 1 \c'{9)\d9 
Js 1 Js 1 
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and LYh = c*((L c * h).ds). Thus we can finally write down the quotient metric 
Gf s ' n (h,k)= f (LY h ,Y k )dx 

JR 2 

(3) = / (L c *h{6), [ F(c(0) - c(0i)) / L c {e l ,6 2 )k{9 2 )ds 2 ds l )ds 

Js 1 x Js 1 Js 1 1 

= [ (L c *h{6),k(6))ds = [[ L c (9,9 1 )(h(8 1 ),k(8))ds 1 ds. 
Js 1 JJs^xs 1 

The dual metric on the associated smooth cotangent space L c * C°°(S 1 ,M. 2 ) is 
similarly: 

Gf S ' n (p,q) = [[ F c (9 7 9 1 )(p(9 1 ),q(9))d Sl ds. 
J Js^s 1 

5.4. The geodesic equation on Emb(S' 1 ,IR 2 ) via conservation of momen- 
tum. A quite convincing but not rigorous derivation of this equation can be given 
using the fact that under a submersion, geodesies on the quotient space are the pro- 
jections of those geodesies on the total space which are horizontal at one and hence 
every point. In our case, the geodesies on Diff(M 2 ) can be characterized by the 
strong conservation of momentum we found above: (p(t)*oj(t) is independent of t. 
If X(t) is the tangent vector to the geodesic, i.e. the velocity X(t) = d t (p o (^^(f), 
then uj(t) is just LX(t) = c*(p(6,t).ds) = c*(p(9,t).\ce(6,t)\.d6) considered as a 
measure valued 1-form instead of a vector field. 

When we pass to the quotient Emb(S' 1 ,K 2 ), a horizontal geodesic of diffeomor- 
phisms tp(t) with ip(0) = identity gives a geodesic path of embeddings c(6,t) = 
ip(t) o 0(0,8). For these geodesic equations, it will be most convenient to take as 
the momentum the 1-form p(0,t) = p(0,t).\ce(0,t)\, the measure factor d9 being 
constant along the flow. We must take the velocity to be the horizontal vector field 
X(t) = F * c(-,t)*(p(6,t).d6). For this to be the velocity of the path of maps c, 
we must have c t (0,t) = X(c(6),t) because the global vector field X must extend 
Ct- To pair p and Ct, we regard p as a 1-form along c (the area factor having been 
replaced by the measure d9 supported on C). The geodesic equation must be the 
differential form of the conservation equation: 



ip(t)*p(-,t) is independent of t. 



More explicitly, if d x stands for differentiating with respect to the spatial coordi- 
nates x, y, then this means: 

dx¥>{t) T \c(e,t)P~(0,t) = enst. 
We differentiate this with respect to t, using the identity: 

d t d x ip(t) = d x (tp t {t)) = d x (X o <p(i)) = (d x (X) o ip(t)) ■ d x ip(t), 

we get 

= d x cp(t) T ■ {{d x (X) T o c{0, t)) ■ p(9, t) + p t (8, t)) . 
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Writing this out and putting the discussion together, we get the following form for 
the geodesic equation on Emb(S' 1 , R 2 ): 

c t (9,t)=X(t)oc(9,t) 

Pt(9,t) = -giadX t (c(9,t),t) -p(9,t) 

X(t) =F*c(-,t) t (p(6,t).d0) 

Note that X is a vector field on the plane: these are not closed equations if we 
restrict X to the curves. The gradient of X requires that we know the normal 
derivative of X to the curves. Alternatively, we may introduce a second vector- 
valued kernel on S 1 depending on c by: 

^(0i,0 2 )=gradF(c(0i)-c(0 2 )). 
Then the geodesic equations may be written: 

ct(6,t) = (F c *p)(9,t) 

where, in the second formula, the dot product is between the two j5's and the vector 
value is given by F' c . 

The problem with this approach is that we need to enlarge the space Diff (M 2 ) to 
include diffeomorphisms which are not C°° along some C°° curve but have a mild 
singularity normal to the curve. Then we would have to develop differential geom- 
etry and the theory of geodesies on this space, etc. It seems more straightforward 
to outline the direct derivation of the above geodesic equation, along the lines used 
above. 



5.5. The geodesic equation on Emb(5 1 ,M 2 ), direct approach. The space of 
invertible pseudo differential operators on a compact manifold is a regular Lie group 
(see [3]), so we can use the usual formula d(A~ x ) = —A~ 1 .dA.A~ 1 for computing 
the derivative of L c with respect to c. Note that we have a simple expression for 
D c mF c , namely 

D c , h F c (6 1 ,6 2 ) = dF{c{6i) - c(9 2 ))(h(9 1 ) - h(9 2 )) = {F' c {9 1 ,9 2 ),h{9 1 ) - h{9 2 )) 
hence 

D c . e L c (9 u 9 2 ) = - f L c (e 1 ,9 3 )D Cth F c (e 3 ,e i )L c (9 4 ,e 2 )d0 3 d9 4 
As 1 ) 2 

= - I L c (9 u 9 3 )((F^9 3 ,9 4 ),£(9 3 ))L c (9 4 ,9 2 )d9 3 d9 4 
Jis 1 ) 2 

+ f L C (9 U 9 4 ) {{F' c (9 4 ,9 3 ),£(0 3 )) L c {9 3 ,9 2 )d9 3 d9 4 

JiS 1 ) 2 

We can now differentiate the metric where 9 = {9\,9 2 ,..., 9 n ) is the variable on 
(S 1 )™: 

D c . e G dm > n (h,k) = f D cJ L c (9 1 ,9 2 )(h(9 2 ),k(9 1 ))d9 

JiS 1 ) 2 



J 



/(Si)= 

L c (9i,9 3 )F' c (9 3 ,9 4 )L c (9 4 , 9 2 ) 
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+ L C (8 U 6 4 )F' C (8 4 , 3 ) L c (6 3 ,d 2 )J(6 3 ))(h(0 2 ), fc(0i)) dd 
We have to write this in the form 

D c . A Gf s ' n {h, k) = Gf s > n (£, H c (h, k)) = G diS > n (K c (t, h),k) 
For H c we use 5(0 5 - 9 3 ) = J L c (8 5 ,8 e )F c (8 e , 6 3 )d8 e = (L c * F c )(8 5 , 8 3 ) as follows: 

D c jG dlS n (h, k) = [ L c (8 5 ,8 6 )(e(8 5 ), ( - L C (8 U 9 3 )F^9 3 ,9 4 )L C (8 4 , 8 2 ) 
+ L c (e 1 ,6 i )F' c (0 A , 8 3 ) L c (8 3 , 2 )) F c (8 6 , 6 3 )(h(6 2 ), fc(*i)>) dd 

Thus 

H c (h,k)(8 ) = f (-L c (e u e 3 )F^e 3 ,e 4 )L c (e 4 ,e 2 ) 

+ L c {8 u 6 4 )F' c (e 4 , 8 3 ) L c (8 3 , fe))*^, W^), d8 

Similarly we get 

D Cti G diS ' n (h,k) = f L c (8 & ,8 h )(F c (8 l ,8 & )(-L c {d u 8 3 )F' c {8 3 ,8 4 )L c (8 4 ,8 2 ) 
Jis 1 ) 6 x 

+ L C {8 U 4 )F' C {8 4 , 8 3 ) L c (8 3 , 6 2 ),£(6 3 ))h(6 2 ),k(6 5 )) dd 

so that 
K c (l,h)(8 ) = 

= f (- (K(9o, 01)^,82), e(6 )) + (F' c {8 , 8 X ) L c (8 1 ,8 2 ),e(8 1 ))) h(8 2 ) d8 
Jis 1 ) 2 v ' 

By 2.4 the geodesic equation is given by 

ctt(Oo) = \H c {ct,ct){e ) - K c {c t ,c t )(6 ) 

Let us rewrite the geodesic equation in terms of L c * c t . We have (suppressing the 
variable t and collecting all terms) 

(L c *ct)t(0o)= / D CtCt L c (e o ,0 1 )ct{e 1 )d0 1 +L c *c u 
Js 1 

= \ J ( F 'c{0iM - K(0o,0i)) {Lc * ct(0 Q ),L c * c t (8\)) Mi 

Since the kernel F is an even function we get the same geodesic equation as above 
for the momentum p(8, t) = L c * Cj = p(8, t).\c$\: 

(1) Pt(0 ) = - [ F^8 ,8 1 )(p(8 ),p(8 1 ))d8 1 

Js 1 

5.6. Existence of geodesies. 

Theorem. Let n > 1. For each k > 2n — | the geodesic equation 5.5 (1) has 
unique local solutions in the Sobolev space of H k -embeddings. The solutions are 
C°° in t and in the initial conditions c(0, . ) and ct(0, . ). The domain of existence 
(int) is uniform in k and thus this also holds in Emb(S' 1 , R 2 ). 
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An even stronger theorem, proving global existence on the level of iJ fe -diffcomor- 
phisms on R 2 , has been proved by [19, 21, 22]. 

Proof. Let c £ H k . We begin by checking that F' c is a pseudo differential operator 
kernel of order — 2n + 2 as we did for F c in 5.3. 

2W R e lia+^^^p+i^n^j^ 

2tt Jr 

where we changed variables as 771 = (c(9i,6 2 ),£) and 772 = (Jc(9i,9 2 ),S,)- So we see 
that F' c {9\,9 2 ) is an elliptic pseudo differential operator kernel of degree — 2n + 2 
(the loss of 1 comes from integrating with respect to 772). We write the geodesic 
equation in the following way: 

c t = F c *q=: Yi(c, q) 

q t = (q,F^q) = J F' c { . ,6){q{9),q{ .)) d9 =:Y 2 (c,q) 

We start with c £ H k where fc > 2n — ^, in the H 2 -open set U k := {c : \cg\ > 
0} C H k . Then q = L c * c t £ H k - 2n+1 and F' c * q £ H k ~ x C H k ~ 2n+1 . By 
the Banach algebra property of the Sobolev space n k - 2n +^ the expression (with 
missuse of notation) Y 2 (c,q) — (q,F[. * q) £ H k ~ 2n+1 . Since the kernel F is not 
smooth only at 0, all appearing pseudo differential operators kernels are C°° off 
the diagonal, thus are smooth mappings in c with values in the space of operators 
between the relevant Sobolev spaces. Let us make this more precise. We claim that 
ci-»Fj*(.)e L(H k , H k+2n ~ 2 ) is C°°. Since the Sobolev spaces are convenient, wc 
can (a) use the smooth uniform boundedness theorem [9], 5.18, so that it suffices to 
check that for each fixed q £ H k the mapping c F' c * q is smooth into ff k + 2n - 2 . 
Moreover, by [9], 2.14 it suffices (b) to check that this is weakly smooth: Using 
the L 2 -duality between H k+2n ~ 2 and H^ k ^ 2n+2 it suffices to check, that for each 
p £ iJ- fe -2«+2 the expression 

J p{9 l )(F' c *q)(9 l )d9 l = 
= J J EM J e i^-02)m 



\c(9 1 ,9 2 )\- 3 .Jv 



-dm) dr h q(9 2 )d9 1 d9 2 



r(i + ttof(I»?iI 2 + M 2 )) 

is a smooth mapping Emb(S' 1 ,R 2 ) — ► R 2 . For that we may assume that c de- 
pends on a further smooth variable s. Convergence of this integral depends on the 
highest order term in the asymptotic expansion in 77, which does not change if we 
differentiate with respect to s. 
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Thus the geodesic equation is the flow equation of a smooth vector field Y = 
(Yi,Y 2 ) on U k x H k ^ 2n+1 . We thus have local existence and uniqueness of the flow 
Fl fe on U k x H k - 2n+1 . 

Now we consider smooth initial conditions Cq = c(0, . ) and q a = q(0, . ) = (L c * 
c t )(0, . ) in C 00 (5' 1 , M 2 ). Suppose the trajectory F\ k (c 0l q ) of Y through these intial 
conditions in U k x H k+1 ^ 2n maximally exists for t £ (—a k ,b k ), and the trajectory 
Fl k+1 (c ,u ) in U k+1 xH k+2 ~ 2n maximally exists for t £ (-a k+1 ,b k+1 ) with b k+1 < 
b k . By uniqueness we have Fl k+1 (c ,u ) = Fl t (co,«o) for t € (— a k +\bk+\)- Wc 
now apply dg to the equation q t = Y 2 (c, q), note that the commutator q [F' c , dg] * 
q = d t h{F' c *q) — F' c * (dgq) is a pseudo differential operator of order — 2n + 2 again, 
and obtain 

de<h = f[F' c ,dg]{.,6){q{e),q{.))de + j F' c {.,6){dgq{6),q{.))d6 
+ J F' c {.,6)(q{6),d e q{.))de 

which is an inhomogeneous linear equation for w = dgq in U k x Jj k + 1 - 2r \ By the 
variation of constant method one sees that the solution w exists in H k for as long 
as (c, q) exists in U x f{ k + 1 ^ 2n j i. e ., for all t £ (~a k ,b k ). By continuity we can 
conclude that w = dgq is the derivative in H k+2 ^ 2n for t = b k +i, and thus the 
domain of definition was not maximal. Iterating this argument we can conclude 
that the solution (c, q) lies in C°° for t £ (—a k , b k ). □ 



5.7. Horizontality for G difi> . The tangent vector h £ T c Emb(5' 1 , M 2 ) is Gf S ' n - 
orthogonal to the Diff (S' 1 )-orbit through c if and only if 

= Gf s > n (h,(x(c)) = [ L c (6 u 6 2 )(h(6 2 ),cg(6 1 ))X(0 1 )d Sl d S2 
Jis 1 ) 2 

for all X £ XiS 1 ). So the G diff '"-normal bundle is given by 

Mf m - n = {h£ CiS 1 ,®?) : (L c *h,v) =0}. 

Working exactly as in section 4, we want to split any tangent vector into vertical 
and horizontal parts as h = h T + h 1 - where h T = X(h).v for X(h) £ X(5' 1 ) and 
where h 1 - is horizontal, (L c * h^-,v) = 0. Then (L c * h,v) = (L c * (X(h)v),v) and 
we are led to consider the following operators: 

Lj,L±:C°°{S 1 )^C°°{S 1 ), 

Lj(f) = (L c *(f.v),v) = (L c *(f.n),n), 

Li(f) = (L c *(f.v),n)=-(L c *(f.n),v). 

The pseudo differential operator Lj is unbounded, selfadjoint and positive on 
L 2 {S 1 ,d9) since we have 

/ Lj(f).fd8= [ (L c (9 1 ,e 2 )f(e 2 )v(e 2 ),f(e 1 ).v(e 1 ))de = \\f.v\\% di!! , n >o. 

Js 1 Jis 1 ) 2 

Thus Lj is injective and by an index argument as in 4.5 the operator Lj is invert- 
ible. Moreover, the operator Lj: is skew-adjoint. To go back and forth between the 
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natural horizontal space of vector fields a.n and the G '"-horizontal vectors, we 
have to find b such that L c * (a.n + b.v) = f.n for some /. But then 

Lj-(a) = -(L c * (a.n),v) = (L c * (b.v),v) = Lj (b) thus b = (Lj)- 1 L± (a). 

Thus a.n+ (Lj)~ 1 L^-(a).v is always G dlff '"-horizontal and is the horizontal projec- 
tion of a.n + b.v for any b. 

Proposition. For any smooth path c in Imm(S' 1 , R 2 ) there exists a smooth path <p 
in Diff(5' 1 ) with <p(0, . ) — Id s i depending smoothly on c such that the path e given 
by e(t, 6) = c(t, ip(t, 6)) is G diS ' n -horizontal: (L c * e u e 6 ) = 0. 

Proof. Let us write e = c o ip for e(t,9) = c(t,ip(t,6)), etc. We look for ip as the 
integral curve of a time dependent vector field 9) on S 1 , given by tpt — £ o ip. 
We want the following expression to vanish. In its computation the cquivariance of 
L c under f G Diff + (S' 1 ) from 5.3(2) will play an important role. 

(L cotp * (d t (coip)),d g (coip)) = (L coip *(c t oip+ ( Cg o <p)<p t ), (c e o <p)<pe) 
= (((L c * c t ) o ip) + ((L c * ((%.£)) ° <P), ( c ° ^Pe) 
= (((L c * c t , c e ) + (L c * (£.c e ), c g )) o ip) ip e . 

Using the time dependent vector field ^ = — (Lj)~ 1 (L c * c t} cg) and its flow ip 
achieves this. □ 

To write the quotient metric on B e , we want to lift normal vector fields a.n 
to a curve C to horizontal vector fields on Emb(S' 1 ,R 2 ). Substituting h = a.n + 
{Ll)- x L^{a).v, k = b.n + (L T )- 1 L ± (b).v in 5.3(3), we get as above: 



The dual metric on the cotangent space is just the restriction of the dual metric 
on Emb(S ,1 ,R 2 ) to the cotangent space to B e and is much simpler. We simply set 
p = f.n, q = g.n and get: 



where x(s) £ R 2 stands for the point in the plane with arc length coordinate s 
and F is the Bessel kernel. Since these are dual inner products, we find that the 
two operators, (a) convolution with the kernel F(x(s) — x(si)).(n(s), n(si)) and (b) 
L~l + Ljr (Lj)^ 1 Ljr are inverses of each other. 

5.8. The geodesic equation on B e via conservation of momentum. The 

simplest way to find the geodesic equation on B e is again to specialize the general 
rule ip(t)*u)(t) = enst. to the horizontal geodesies. Now horizontal in the present 
context, that is for B e , requires more of the momentum u)(t). As well as being 
given by c*(p(s).ds), we require the 1-form p to kill the tangent vectors v to the 
curve. If we identify 1-forms and vectors using the Euclidean metric, then we may 
say simply p(s) = a(s).n, where a is a scalar function on C. But note that if you 





F(x(s) - a;(si)).(n(s), n(s 1 )).f(s)g(i).dsds 1 
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take the momentum as c*(a(s)n(s)ds) and integrate it against a vector field X, 
then you find: 



where ix is the 'interior product' or contraction with X taking a 2-form to a 1-form. 
Noting that 1-forms can be integrated along curves without using any metric, we 
see that the 2-form along c defined by {a(s).(dx A dy) c r s \} can be naturally paired 
with vector fields so it defines a canonical measure valued 1-form. Therefore, the 
momentum for horizontal geodesies can be identified with this 2-form. 

If <p(x, t) is a horizontal geodesic in Diff (K 2 ), then the curves Ct = image(c(-, t)) 
are given by C t — ip(Co,t) and the momentum is given by a(0,t).(dx A dy), where 
c(6, t) parametrizes the curves C t - Note that in order to differentiate a with respect 
to t, we need to assign parameters on the curves C t simultaneously. We do this in 
the same way we did for almost local metrics: assume eg is a multiple of the normal 
vector nc- But 6 ip(9 ,t) gives a second map from Co to C t : in terms of 9, 
assume this is 6 = <fi(6o, t). Then the conservation of momentum means simply: 



Let X be the global vector field giving this geodesic, so that tpt = X o ip. Note that 
(ft = {{X oc, v) /\cg\)o(p. Using this fact, we can differentiate the displayed identity. 
Recalling the definition of the flow from its momentum and the identifying TcB e 
with normal vector fields along C, we get the full equations for the geodesic: 



Note, as in the geodesic equations in 5.4, that we must use F to extend X to the 
whole plane. In this case, we only need (a) the normal component of X along C, 
(b) its tangential component along C and (c) the divergence of X along C. These 
are obtained by convolving a{s) with the kernels (which we give now in terms of 
arc-length) : 




a(ip(6o,t),t).det(D x <p)(c(6o,0),t) is independent of t. 



X 



(X, n) ■ n 

-(X, v)D s (a)-div(X).a 
F * c*(a(s)n(s)ds) 



*r(*i>*2) 

F v c n {s llS2 ) 



F(c(si) - c(s 2 ))(n(si),n(s 2 ))ds2 
F(c(si) - c(s 2 ))(v(s 1 ),n(s 2 ))ds 2 
(gradF(c(si) - c(s 2 )),n(s 1 ))ds 2 



Then the geodesic equations become: 



C t = (F c n ™ * a).n 

a t = -(F c OTl * a)D s (a) - (F c div * a). a 



Alternately, we may specialize the geodesic equation in 5.4 to horizontal paths. 
Then the v part vanishes identically and the n part gives the last equation above. 
We omit this calculation. 
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6. Examples 

6.1. The Geodesic of concentric circles. All the metrics that we have studied 
are invariant under the motion group, thus the 1-dimensional submanifold of B e 
consisting of all concentric circles centered at the origin is the fixed point set of the 
group of rotations around the center. Therefore it is a geodesic in all our metrics. 
It is given by the map c(t, 6) = r(t)e l6 . Then eg = ire 10 , v c — ie 10 , n c = —e 10 , 
1(c) = 2nr(t), n c = i and c t = r t e ie = -r t .n c . 

The parametrization r(t) can be determined by requiring constant speed a, i.e. 
if the metric is G(h, k), then we require G c (ct, Ct) = r 2 G c (n c , n c ) — a 2 , which leads 
to y/ G c (ct, ct) dt — ±yjG c (n c , n c ) dr. To determine when the geodesic is complete 
as r — > and r — > oo, we merely need to look at its length which is given by: 



\J G c (c t ,c t )dt = / \J G c (n c , n c )dr, 
o Jo 

and we need to ask whether this integral converges or diverges at its two limits. 

Let's consider this case by case. 

The metric G*: The geodesic is determined by the equation: 

G*(ct, ct) = 2nr ■ $ ^7rr(i), ■ r\ = a 2 . 

Differentiating this with respect to t leads to the geodesic equation in the standard 
form r t t — r 2 f(r). It is easily checked that all three invariant momentum mappings 
vanish: the reparameterization, linear and angular momentum. 

Theorem. If $(27rr, 1/r) w C\r a (resp. C2r b ) as r — > (resp. oo), then the ge- 
odesic of concentric circles is complete for r — > if and only a < — 3 and is 
complete for r — > oo if and only if b > —3. In particular, for ip = £ k , we find 
k = a = b and the geodesic is given by r(t) = cnst.t 2 ^ k+3 . For the scale invariant 
case $(t, k) — ^p- + we find a = b = —3, the geodesic is given by r(t) — e^ 2 " 1 
and is complete. Moreover, in this case, the scaling momentum is constant in 
t along the geodesic. 

The proof is straightforward. 

The metric G imm - n Recall from 4.1 the operator L n , c = I + (-l) n A.D 2n . For 
c(t,9) = r(t)c w 6.1 we have 

LM = (1 + (-l) n ^)(r t e id ) = r t (l + A) c " 

which is still normal to c$. So f ^ c(t, . ) is a horizontal path for any choice of 
r(t). Thus its speed is the square root of: 

/ a \ 

™m,n/„ „ \ _ o_„ I 1 i I „2 2 



G ~ {ctiCt ) = 2vr-[l + —j-ri = o- 

For n = 1 this is the same as the identity for the metric with n) = 1 + An 2 
which was computed in [12], 5.1. An explanation of this phenomenon is in [12], 3.2. 

Theorem. The geodesic of concentric dries is complete in the G lmm '™ metric if 
n > 2. For n = 1, it is incomplete as r — > but complete if r — > oo. 
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The metric G d To evaluate the norm of a path of concentric circles, we now 
need to find the vector field X on R 2 gotten by convolving the Bessel kernel with 
the unit normal vector field along a circle. Using circular symmetry, we find that: 



X(x,y) = f(r)(-J)r t 
\r r I 



I — A(d rr + -d r )\ f = except on the circle r = r 



r 



f e C everywhere, f(jo) = 1 



For n — 1, we can solve this and the result is the vector field on R 2 given by the 
Bessel functions h and K\: 

Ml/vg) if r < rn 

/i(ro/V3) " 

Using the fact that the Wronskian of h,K\ is 1/r, we find: 

G dis -\r t n,r t n) = [ ((I - AA)X, X)r\ ^ 



X(x,y) 



Using the asymptotic laws for Bessel functions, one finds that the geodesic of con- 
centric circles has finite length to r = but infinite length to r = oo. 

For n > 1, it gets harder to solve for X. But lower bounds are not hard: 
G diff <"(n, n) = inf f ((I - AA) n X, X) 

X,(X,n) = l onC r J 

>A n inf f(A n (X),X) = M(r) 

/ X,(X,n) = l on C r \ J dcf 
yX^O, when x^oo J 

Then M(r) scales with r: M(r) = Af(l)/r 2 "~ 2 , hence the length of the path when 
the radius shrinks to is bounded below by ^dr/r 11 " 1 which is infinite if n > 1. 
On the other hand, the metric G dlS ' n dominates the metric G fdlff ' 1 so the length of 
the path when the radius grows to infinity is always infinite. Thus: 

Theorem. The geodesic of concentric dries is complete in the G dlii ' n metric if 
n > 2. For n = 1, it is incomplete as r — > but complete if r — > oo. 

6.2. Unit balls in five metrics at a 'cigar'-like shape. It is useful to get a sense 
of how our various metrics differ. One way to do this is to take one simple shape 
C and examine the unit balls in the tangent space TcB e for various metrics. All of 
our metrics (except the simple L 2 metric) involve a constant A whose dimension is 
length-squared. We take as our base shape C a strip of length L, where L 3> V~A, 
and width w, where w -C y/A. We round the two ends with semi-circles, as shown 
in on the top in figure 1. 

As functions of a normal vector field a.n along C, the metrics we want to compare 
are: 

(1) G£(a,a) = J c (l + AK 2 )a 2 .ds 

(2) G™ m!l (a, a) = inf J c (|a • n + b ■ v)\ 2 + A\D s (a -n + b- v)\ 2 ^ds 
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C 



2w ^ involution ' 



c 



Iff 



il± Cji , ' 



c 



3 



Figure 1. The cigarlike shape and their deformations 



(3) G d f-\a,a) = ^ inf JJ R2 (|X| 2 + A\DX\A dx dy, 

^ {X,n)=a ) 

(4) G ^ 2 (a,a) = ^ inf JJ R2 (|X| 2 + 2A\DX\* + A^X 

(X,n)=ct ) 



\ 2 ^j dxdy 



The term in the last 2 metrics is put there so that the double integrals have 
the same 'dimension' as the single integrals. In this way, all the metrics will be 
comparable. 

To compare the 4 metrics, we don't take all normal vector fields a.n along C. 
Note that C has an involution, which flips the top and bottom edges. Thus we 
have even normal vector fields and odd normal vector fields. Examples are shown 
in figure 1. We will consider two even and two odd normal vector field, described 



below, and normalize each of them so that J c 



1. They are also shown in 



figure 1. They involve some interval / along the long axis of the shape of length 
A>id. The interval determines a part I t of the top part of C and lb of the bottom. 



(1) Let a = +I/V2A along I t and a = — 1/V2A along h, a zero elsewhere 
except we smooth it at the endpoints of I. Call this odd vector field if- . 

(2) Fix a high frequency / and, on the same intervals, let a(x) = ± sin(/.a;) 
Call this odd vector field ipf. 

(3) The third vector field is even and is defined by a(x) = \Jl^{n, at the 
right end of the curve, being zero along top, bottom and left end. Call this 
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tp x - The factor in front normalizes this vector field so that its L 2 norm is 
1. 

(4) Finally, we define another even vector field by a — +1/V2A on both I t and 
lb , zero elsewhere except for being smoothed at the ends of I. Call this ip+ . 

The following table shows the approximate leading term in the norm of each of 
these normal vector fields a in each of the metrics. By approximate, we mean the 
exact norm is bounded above and below by the entry times a constant depending 
only on C and by leading term, we mean we ignore terms in the small ratios w/y/A, 
X/VA: 

function G A G™" 1 ' 1 G^' 1 G diff < 2 



<Pf 

fx 



1 
1 
1 

A/w 2 



1 

(VAf) 2 
1 

VA/w 



1 

VAf 
y/A/w 

w \og(A/w 2 ) 



l 



(Va/) 3 

{yfA/wf 
y/A/w 



Thus, for instance: 



G A (lfix,(Px) = — 



{ l + A ^ ){n oy ds 



right end 

.-2> 



= ^(1 + Aw- 2 )£(right cnd)Avc((n, £) 2 ) 
= (1 + Aw- 2 ) w A/w 2 

The values of all the other entries under G are clear because k = in their support. 

To estimate the other entries, we need to estimate the horizontal lift, i.e., the 
functions b or v. To estimate the norms for G imm . 1 j wc take b = in all cases except 
(fx and then get 

Gimm.l/ \ II n2 
C ( a ' a ) = \\ a \\H\ 

the first Sobolev norm. For a — iff, we simplify this, replacing the full norm by 
the leading term A(D s (a)) 2 and working this out. To compute G lmm -' 1 (ip x , ( Px), let 

k = \J ^ be the normalizing factor and lift a.n along the right end of C to the M 2 
vector field k.-^. This adds a tangential component which we taper to zero on the 
top and bottom of C like fc.e -1 /^. This gives the estimate in the table. 

Finally, consider the 2 metrics G ' k , k — 1,2. For these, we need to lift the 
normal vector fields along C to vector fields on all of ffi. 2 . For the two odd vector 
fields / = <p- and / = <pf, we take v to be constant along the small vertical lines 
inside C and zero in the extended strip —w<y<w,x^I and we define v outside 

— w < y < w by: 



v(x,-y-w) = F(x,y)-§jj, 

A:Vl(l + ^ 2 ) fc - 1/2 ./(0 
ir{l + A({ 2 + V 2 )) k 



v(x, y + w) 

We check the following: 

(a) ((7 - AA) k F) A = + M 2 ) k ~ 1/2 f(Z) is inde P- of V hence 

support(7 - AA) k F) c {y = 0} 
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(b) J Fdrj = /, hence F\ y=0 = f. 

Thus: 

/) « ^ /j[ 2 ((/ - AA) fc ^, F)dxdy 

= I /a + ^ 2 )*- 1/2 i/(oi 2 de - f ii/ii^-vs 

The leading term in the fc 4/l Sobolcv norm of iff is (v A A/) 2fe , which gives these 
entries in the table. 

To estimate Gq '^{if+, we define w by extending <p + linearly across the 
vertical line segments — w < y < w,x € I, i.e. to ip + (x)y/w. This gives the leading 
term now, as the derivative there is ip+/w. In fact for any odd vector field a of 
L 2 -norm 1 and for which the derivatives are sufficiently small compared to w, the 
norm has the same leading term: 

Gc s -\a,a) w VA/w. 

To estimate G^ S ' 2 (ip + , tp+), we need a smoother extension across the interior of 
C. We can take <p + (x).|(^ — |(^) 3 )- Computing the square integral of the second 
derivative, we get the table entry G^ S ' 2 (a,a) w (y/A/w) 3 . 

To estimate G^ S ' k+1 (<p x , ip x ), we now take v to be 



c(k, A, w) 



d_ 

dx 



where D is the disk of radius w containing the arc making up the right hand 
end of C, and where c(k,A,w) is a constant to be specified later. The function 

— (M^j K k {^L} is the fundamental solution of (I - AA) k+1 and is C 1 for 
k > but with a log pole at for k = 0. Thus: 

(7 - AA) k+1 v = 2irk\Ac(k, A, w)xd--^ 

while, up to upper and lower bounds depending only on k, the restriction of v to the 
disk D itself is equal to log(vC4/w)c(0, A, w)w 2 if k = and simply c(k, A, w)w 2 for 
k > 0. By symmetry v is also constant on the boundary of D and thus v extends 
ip x if we take c(0,A,w) = c / \og(y^A/w)w 5 ^ 2 if k = and c(fc, A, to) = c^/w 5 / 2 if 
fc > (constants Ck depending only on k). Computing the 7f fc -norm of v, we get 
the last table entries. 

Summarizing, we can say that the large norm of ip x is what characterizes G A ; 
the large norms of <p + characterize G dlff ; and the rate of growth in frequency of the 
norm of iff distinguishes all 4 norms. 
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